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Abstract 

In  this  paper  we  consider  initial-boundary  value  problems  for 
an  energy  conserving  dispersive  hyperbolic  equation,  the  KLein- 
Gordon  equation.   This  equation  contains  the  main  feature  of  dispersion: 
The  speed  of  propagation  depends  on  the  frequency. 

The  primary  purpose  of  this  paper  is  to  compare  the  asymptotic 
expansion  of  solutions  obtained  by  a  technique  which  we  call  the  ray 
method  with  the  asymptotic  expansion  of  the  exact  solution.   In  every 
case  considered,  the  solutions  agree. 

Solutions  are  obtained  for  a  series  of  initial-boundary  value 
problems  in  one  space  dimension  with  variable  coefficients.   As  a 
by-product  of  our  work  we  find  a  new  feature  which  we  call  space -time 
diffraction.   This  phenomenon  has  the  following  physical  interpretation: 
A  portion  of  the  energy  of  a  wave  reaches  a  boundary  surface  and  then 
gradually  leaks  off,  leaving  a  diminishing  residue  on  the  boundary 
for  all  time. 


Introduction 

In  the  past  few  years  there  has  been  a  great  deal  of  Interest  in  the 
asymptotic  expansion  of  solutions  of  boundary  value  problems  and  initial- 
boundary  value  problems  for  partial  differential  equations.   These  expansions 
are  obtained  for  large  distance  and/or  time  or  for  large  values  of  a  parameter 
which  may  appear  "naturally"  in  the  equation  and/or  initial  conditions,  boundary 
conditions,  etc. 

Many   of  these  problems  can  be  treated  by  a  method  in  which  special  curves, 
called  "rays",  in  space  or  space-time  plan  an  essential  role.   The  asymptotic 
solution  is  determined  in  terms  of  "amplitude"  and  "phase"  functions  which  satis- 
fy ordinary  differential  equations  on  these  rays.   These  equations  can  often  be 
solved  explicitly,  so  that  the  solution  to  the  original  problem  is  known  if  tiie 
"initial  -values"  on  the  rays  of  the  amplitude  and  phase  functions  are  known.   In 
some  cases  these  initial  values  can  be  determined  directly  from  the  data  (i.e., 
initial  conditions,  boundary  conditions,  etc.)  of  the  given  problem,  while  in 
other  situations  an  indirect  method  is  necessary.   This  method  involves  the 
exact  solution  of  a  "canonical  problem",  i.e.,  a  problem  simpler  than  that  given 
originally  but  having  the  same  local  features  as  the  given  problem. 

This  "ray  method"  has  been  studied  extensively  by  J.B.  Keller  and  his 
associates  at  New  York  University  for  the  reduced  wave  equation  and  related 
elliptic  differential  equations.  A   significant  advantage  of  the  method  is  that 
it  eliminates  the  necessity  of  first  finding  the  exact  solution  in  order  to  find 
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the  asymptotic  expansion.  Consequently,  many  problems  for  which  the  exact 
solution  is  not  known,  can  be  treated  by  the  ray  method.   In  particular,  it 
has  led  to  a  "geometrical  theory  of  diffraction"  [  1  ]. 

Recently,  the  ray  method  has  been  extended  to  time  dependent  problems 
and,  in  particular,  to  dispersive  hyperbolic  equations  [2],   One  example  of 
this  class  is  the  Klein-Gordon  equation  which  may  be  used  to  describe  the 
propagation  of  electromagnetic  waves  in  certain  plasmas  [2].   Initial-boundary 
value  problems  for  this  equation  are  presently  being  studied,  since  it  is  the 
simplest  example  of  a  large  class  of  dispersive  hyperbolic  equations  which 
conserve  energy.   It  clearly  exhibits  the  main  feature  of  dispersion:   The 
speed  of  propagation  varies  with  frequency. 

As  of  now  there  is  no  general  proof  that  the  ray  method  yields  the 
asymptotic  expansion  of  the  exact  solution.   Consequently,  a  major  goal  of 
the  present  research  is  to  "check"  the  solutions  obtained  by  this  method  by 
comparing  them  with  the  asymptotic  expansion  of  the  exact  solution  in  special 
cases  where  the  latter  can  be  obtained.   Such  checks  have  previously  been  ob- 
tained for  the  reduced  wave  equation  [  1,3,U  ]•   Here  they  are  undertaken 
for  initial-boundary  value  problems  for  the  KLein-Gordon  Equation.   In  every 
case  examined  the  asymptotic  solution  obtained  by  the  ray  method  agrees  ex- 
actly with  the  expansion  of  the  exact  solution. 

¥e   shall  consider  the  Klein-Gordon  equation  in  one  space  dimension 
with  non-constant  coefficients.   Initial-boundary  value  problems  for  the 
homogeneous  equation  are  solved  by  the  ray  method  and  the  asymptotic  expan- 
sion of  the  Green's  function  is  also  obtained.   The  distinguishing  feature 
in  the  initial-boundary  value  problems  considered,  is  the  nature  of  the  de- 
pendence of  the  initial  data  on  the  large  parameter. 


A   method  for  obtaining  the  asymptotic  expansion  of  the  exact  solu- 
tion is  also  developed  and  applied  to  each  of  the  problems  for  which  a 
solution  was  obtained  by  the  ray  method.  As  previously  noted,  the  solu- 
tions agree  in  each  case. 

As  a  by-product  of  our  work  we  find  a  new  feature  which  we  call 
space-time  diffraction.   This  phenomenon  occurs  in  a  region  of  space-time 
in  which  the  leading  term  of  the  asymptotic  expansion  is  zero.   The  name, 
space-time  diffraction,  is  suggested  by  the  similarity  to  an  analogous 
problem  for  the  reduced  wave  equation  [  3J  tut  it  differs  in  that  it  is  a 
time-dependent  phenomenon,  with  the  following  physical  interpretation:  A 
portion  of  the  energy  of  a  wave  reaches  a  boundary  surface  and  then  gradually 
leaks  off,  leaving  a  diminishing  residue  on  the  boundary  for  all  time. 


(l)  Introduction  to  the  Ray  Method 

The  partial  differential  equation  to  be  studied  is 

(1)     u  -  u^^  -  X  b  (x)u  =  0. 
^         XX    tt 

We  seek  the  first  term  of  the  asymptotic  expansion  of  u(t,x)  for  large  \. 

(l.l)  The  Dispersion  Equation  and  the  Transport  Equation 
The  solution  will  "be  assumed  to  be  of  the  form 


(2)  u(t,x)  =  exp 


iXs(t,x)   z(t,x)l  fl  +  0(^)1 


When   (2)   is   substituted  into   (l)  and  the  coefficients  of  X     and  X  respectively 
are   set   equal  to  zero,   the  resulting  equations  are 

(3)        s^-    s^   +b2(x)   =  0 


and 

(k)      2{s jz.^-    s^z^)    +  zD    s   =  0,      □    =         —5 


XX         t    f 


^  at 


We  define 


(5)        k   =  s^,    u)   =  -s^ 


and  then   (3)  becomes 


(6)       k^-  co^   +  b^(x)   =  0. 


If  b(x)  were  constant  and  the  plane  wave  solution  to  (l)  were  sought, 
one  would ' obtain  (6)  immediately.   This  algebraic  equation  is  called 
the  dispersion  relation.   The  first  order  partial  differential  equation,  ■ 
(3)  (which  is  equivalent  to  (6)  ),  will  be  called  the  dispersion  equation. 
(k)   will  be  called  the  transport  equation. 

(1.2)  The  Solution  of  the  Dispersion  Equation 

The  method  of  characteristics  [p]]  will  be  employed  to  solve  the 
dispersion  equation.   The  characteristic  equations  are 

(7)  X  =  2k,  t  =  2co,  _w  =  0,  k  =  -2b(x)b'(x). 

The  solution  of  this  system  of  ordinary  differential  equations  defines  curves 
in  (t,x)  space,  which  will  be  called  rays.   The  latter  pair  of  equations 
in  (7)  are  used  to  determine  w  and  k  on  the  rays.   The  equation  for  u 
yields 

(8)  w  =  constant. 

Equation  (6)  then  gives  k  =  k(x,w).  With  x  as  independent  variable  on  the 
rays,  the  time,  t,  and  the  phase,  3,  can  be  shown  to  satisfy  the  following 
ordinary  differential  equations  by  using  (5)  and  (7): 


(q)  ^  =       .^ 

C^n^         d.s  _     ,      dt    .    ,  dt   * 
dx    X     t  dx         dx 

*  If  (9)  and  (6)  were  used  in  (10)  we  would  have 

dx   i^^^^^^'^  k(3^;^  -  k(3^7J)      -  -  k(x,co)   • 

ds 
This  result  has  the  virtue  of  expressing  -=—  completely  in  terms  of  the  in- 


dependent variable  on  the  ray.   However,  for  comparison  later  in  the  paper, 
the  form  of  equation  (lO)  best  serves  our  purpose. 


The  expression  ^    ' — -,   obtained  by  inverting  (9),  will  be  referred  to  as 
the  group  speed. 

(1.2.1)  Ray  Construction  and  Phase  for  a  One  Dimensional  Initial  Ivfe.nifold. 

If  s  were  given  on  some  curve  in  (t,x)  space,  one  could  in  general 
find  s(t,x)  for  all  (t^x).   As  a  special  case,  let  us  suppose  that 
s(0,x)  =  s  (x)  is  given.   This  initial  data  for  s  could  be  described 
parametrically  by 

(11)  t  =  0,  X  =  a;   s  =  s^(cT). 

The  method  used  to  obtain  s   (0)   from  a  given  problem  for  u(t,x)  will  be 
explained  later  in  the  discussion.   Here  we  shall  only  show  how  s(t,x) 
is  obtained  once  s  (o)   is  known.   From  (5)  and  (ll) 

(12)  k(a,oj)  =  s^(a). 
Then  (6),    (8)  and  (12)  yield 


(13)   0)  =  h(s^(a),  0)  =  t  hjs'ja),    a);     h^(k,a)  =  /k^  +  b^(a) 


*     We  note  at  this  point,  the  following  convention.   When  a  specific  choice 

of  sign  is  dictated,   -th  or  -h  will  be  written,  but  w^  will  be  written 

instead  of  h2  . 
o 
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With  cj  given  lay  (13)  k  can  be  determined  from  (6)  up  to  a  choice  of 
sign.  By  requiring  (12)  to  "be  satisfied  at  x  =  a  this  ambigiiity  is 
eliminated  yielding  the  solution 

uj  -  b  (x). 
The  solution  to  (9)  for  which  t  >  0  is  given  by 

nS         h^ 

(15)       *  "  )     /  2°  2  ^  '^'^'   \  "  min(x,a),  x  =  max(x,a), 
^x^  /cj  -  b  (ti) 


With  this  result  and    (lU)  the   solution  to   (10 )  becomes 

r  ^ 

(l6)  s(t,x)    =   s   (a)    +   sgn  to  <  \       /u   -  b   (i])      dr]   -  h  t 


^< 


To  obtain  (15)  and  (l6)  we  have  also  used  the  fact  that  sgn(x-a)sgn  k 

+ 
=  sgn  00  which  is  a  consequence  of  the  requirement  that  t  >  0.   If  s"  are 

the  solutions  associated  with  sgn  oj  =  -  1  then  (l6)  becomes 
(17)       s*(t,x)  =  sJa)±M  /J^^U)      dTi-h^tV. 

Equation  (15)  yields  a  one  parameter  (cr)  family  of  rays  and  (17)  gives 

2   2 
s  in  terms  of  the  parameter.   (Note  that  w  =  h  is  given  by  (13).  ) 

(1.2.2)    Ray  Construction  and  Phase  for  a  Zero  Dimensional  Initial  Manifold 

Another  problem  for  which  the  solution,  s(t,x)j  will  be  required  is 

the  problem  where  s  is  given  only  at  a  point  (0,x  ).   The  solutions  of 

interest  will  be  those  for  which  the  rays  emanate  from  the  point  (0,x  ). 

We  assume  the  given  data  to  be  s(0,x  )  =  s  . 


From  the  dispersion  relation  it  follows  that 


(Id) 


k(x  ,uj)  =  +  /7~?0O  ;  k(x,w)  =  +  fF~^{x) 


The  signs  of  k(x  ,^)   and  k(x,w)  are  always  chosen  so  that  they  agree  for 
any  particular  ray.   If  a   is  replaced  "by  x^  and  s^(a)  is  replaced  by 
s  in  (15)  and  (17)       the  equations  thus  obtained  will  be  the 

equations  of  the  rays  and  phase  for  the  zero  dimensional  initial  manifold. 

2    2 
In  this  case  w  is  a  free  parameter,  "  >b  (x^).   As  w  varies,  the  initial 

value  of  the  group  speed  varies  from  -1  to  1.   Note  that  the  comment 

following  (17)  is  no  longer  valid  —  a  is  not  the  parameter  on  the  rays  — 

and  the  initial  value  of  k  is  now  given  by  (15). 

(1.3)   Solution  of  the  Transport  Equation 

In  [2J  it  is  shown  that  the  transport  equation  can  be  written  as  an 
ordinary  differential  equation  along  the  rays.   To  apply  the  results  of 
that  paper  to  the  discussion  below  it  is  assumed  that  fl5);  the  equation 
for  the  rays,  is  invertible  and  gives  x  as  a  function  of  t  and  the  parameter. 
It  is  also  shown  in  [2]  that  z  at  time  t  is  related  to  z  at  time  t.,  by 


(19)         z(t)  =  z(t^) 


*     These  equations,  with  the  indicated  substitutions,  are  given  by  (1.30)  and 
(1.31)  below. 


If  the  rays  are  viewed  as  equations  defining  a  transformation  of  points 
of  a  parameter  space  to  points  x(t)  for  fixed  t,   then  j(t)  is  the  jacobian 
of  that  transformation: 


(20) 


j(t) 


^ 
^ 


or  j(t)  = 


^ 
^ 


The  quotient  — .  /  .,   is  the  jacobian  of  transformation  by  rays  of  points 
X  at  time  t,  to  points  x  at  time  t. 

(1.3.1)  The  One  Dimensional  Initial  Ifenifold 

It  is  assumed  that  for  t  =  0,  x  =  a;  z(0,a)    =  z    (a).      Again,  the 
explanation  of  how  the  initial  value  is  obtained  will  be  postponed.   In 
(15),  cr  is  taken  as  parameter  and  w  =  _  h  (a).   The  equation  is  differentiated 

implicitly  with  respect  to  a  with  t  fixed,  and  one  obtains  a  linear  equation 

Sx 

in  -^  .   We  note  that  there  are  three  terms  in  this  equation  since  the 

upper  and  lower  limit  of  integration  and  the  integrand  are  functions  of 
CT.   This  equation  yields 


(21)   j(t(x))  = 


ax 


OJ 


k3(Tl,a)) 


<3.r\ 


The  jacobian,  j(t(x))  ,  as  given  by  (21;,  might  be  zero  for  some  value  of  x. 
There  might  even  be  a  curve  in  (t,x)  space  on  which  j(t(x))  =0.   Such  curves 
or  Points  are  called  caustics.   Their  effect  on  the  solutions  is  discussed  in 
[2J .   For  simplicity  it  is  assumed  that  there  are  no  caustics.   If  there  were 
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caustics,  our  results  would  only  be  valid  for  times  prior  to  the  arrival 
time  at  the  caustic. 

For  (t,x)  =  {0,a)  ,    J  =  1  and  k  =  s^(a).   In  (19)  we  take  t  =  0, 
t(x)  given  by  (l5)  and  j(t(x))  given  by  (21).   This  yields 


(22)       z(t(x))  =  z(o) 


J(t(x)) 


(1.3' 2)   The  Zero  Dimensional  Initial  Manifold 

From  section  (1.2.2),  to  must  now  be  the  parameter.   In  this  case, 
the  rays  emanate  from  a  point  (0,x  ),  so  that  the  jacobian  of  transformation 
must  be  zero  there.   If  z(t(x))  were  finite  at  t  =  0,  setting  t  =  0  in  (19) 
shows  that  z(t(x))  would  be  identically  zero  for  t  >  0.   In  order  that  z(t(x)) 
be  finite  and  non-zero  for  t  >  0,  it  must  have  a  singularity  at  t  =  0  of  one 
half  the  order  of  the  zero  of  j(t).   This  order  of  vanishing  can  be  determined 

directly  from  the  ray  equations.   We  assume  that  j  can  be  expanded  in  a  Taylor 

Sx 
series  with  constant  term  zero.   Then  if  ^r—  >  0  for  1  >  >  t  >  0  we  may  find 

OU) 

the  leading  term  by  noting  that 

(23)      ^U(t)]=f(i^).^  (^). 

—  is  given  by  (9)  with  k  given  by  (lU).   We  may  carry  out  the  indicated 
operations  in  (23)  and  evaluate  at  (0,x  ),  obtaining 
(21+)       j(t)  =  j^t  +  O(t^); 

(25)       J^  =  |b  (x^)w"  k"  (x^,w)| 

Therefore,  for  the  solution  to  be  finite  and  non-zero,  z  must  be  of  order 

t  ^  in  the  neighborhood  of  t  =  0.   We  define 


(26)       z  =  lim  z(t)t2 
t  ^  0 

and  then,  using  (19)^ 
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C^)       z(t(x))  =z  [-3^]^ 


In  order  to  find  j(t(x))  in  this  case,  the  procedure  of  section 
(1.3.1)  is  repeated,  but  with  a   replaced  by  x  and  k(x,u))  defined  by  (I8); 
(15)  is  differentiated  implicitly  with  respect  to  w  with  t  fixed.   Two  terms 
are  obtained— the  upper  limit  and  the  integrand  of  (15)  are  fimctions  of  co. 
The  solution  for  j(t(x))  obtained  by  this  method  is 


(28)       j(t(x))  =  f^ 


'LO 


o      ^x^   (a)^-b^(Ti)) 


ax  _/u/     -'-'-^    ^^  ^2 


£Mr__b:(jii 


^n        K'^~TJ~XJ^  ^'^' 


Equations  (25)  -  (28)  allow  us  to  express  z  as  a  fixnction  of  x  and  u) 
on  the  rays. 
(1.4)  An  Example 

As  an  example,  the  leading  term  of  the  asymptotic  expansion  of  the 
Riemann  function  for  equation  (l)  will  be  found.   The  initial  data  for 
this  problem  is 


(29)        u(0,x^)  =  0  u^(0,x^)  =  5(x-x^), 


3sume 


Since  the  initial  data  for  u(t,x)  vanishes  except  for  x  =  x  ,  we  as; 

it 

that  all  our  rays  emanate  from  the  point  (0,x  ).   Therefore,  the  results 
of  sections  (1.2.2)  and  (1.3.2)  will  be  used. 

The  rays  and  phase  respectively  are  given  by 


(30)   "t  =  /  -  dT]  ;  X  =  min(x,x  ),  x,^  =  max(x,x  ) 

4  /XT^ 


and 


^  W  r  ^  .rr—^, 


(31)   s-(t,x)  =  s^  ^  \  /  J"^  -   ^  ("H)  ^^  -  \^ 


*  This  assumption  shall  be  justified  when  the  solution  obtained  is  compared  to 
the  asymptotic  expansion  of  the  exact  solution,  obtained  in  section  (2). 
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These  formulas  are  obtained  from  (15)  and  (17)  by  replacing  o  by  x  and 
s  (a)  by  s~.      The  choice  of  sign  in  (3I)  is  not  determined  for  the  zero 
dimensional  initial  manifold.   Since  the  equation  is  linear,  the  solution 
can  be  obtained  by  superposition: 

(32)         u(t,x)~  ^  z-   (t,x)  exp[iAs*(t,x)] 

+ 

with  s~  given  by  (3I)  and  z~  given  by  (27)  with  z   replaced  by  z~. 

If  b  were  constant,  we  could  carry  out  all  indicated  integrations  in 
(28),  (30)  and  (31)  to  obtain  the  solution 


(33)        u~ 


r2*t-l/2exp[ix[s^l    (/ 


'ca)  -b    |x-x  I  -h  t 
'    o'   o  ; 


^o  /j2_^2 

t  =  Ix-x  I   ;     Ix-x  I  =  -^^ t. 

o'   r-r — -   '     I    o'     h 


r2  ,  2  '^  "o 


'u)  -b 


u(t,x)  is  given  parametrically  with  parameter  oj  by  (33)  except  for 

+      + 

the  constants  z~  and  s~.   For  constant  b  the  exact  solution  can  also  be 

o 

obtained.   If  this  solution  is  expanded  asymptotically  and  the  expansion  is 
compared  with  (33  )j  z~  and  s~  can  be  detennined.  We  now  assume  that  these 
values  are  the  same  for  the  case  of  variable  coefficients  if  the  constant 
b  is  replaced  by  b(x  ).   R.  M.  Lewis  calls  this  technique  for  obtaining 
z  and  s   the  indirect  method.   The  asymptotic  expansion  of  the  exact 


om 


solution  for  constant  coefficients  has  been  carried  out  in  [2]  •   Fr 
equations  (62)  and  (63),  section  (3)  of  that  paper. 


(3h)  s^    =   0; 

The  solution  for  variable  b,  obtained  by  using  (27),  (3I),  (32)  and  (3^+), 
is  given  by 
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X 

> 


{4f     /"^-     l3^(Tl)dTl    -     \t     +^ 


cos^x|    /     /u  -  b    (ii)dTi  -   h^t 


< 


^\ 


(35)  ^-r    , .  s — ^ —       >/2 


^< 


(cj  -b    (ti))- 


V 


Equations  (30)  and  (35)  give  the  asymptotic  expansion  of  u(t,x) 
parametrically  with  parameter  w. 

In  section  (2)  the  asymptotic  expansion  of  the  Riemann  Function  will 
be  obtained  from  the  exact  solution.   The  assumptions  made  to  obtain  the 
ray  method  expansion  will  be  justified  by  the  fact  that  the  solutions  agree. 

(1.5)  The  Turned  Rays  and  the  Field  Associated  With  Them 

We  assume  that  b(x)  is  a  monotonically  decreasing  function:  b'(x)<  0. 
Then  for  each  value  of  OJ  there  exists  at  most  one  point,  x  ,  for  which 

(jl) 

2         2 
b  (x  )  =  w  .   A  point  of  this  type  will  be  called  a  turning  point.   In 

this  section  the  ray  method  will  be  studied  in  the  neighborhood  of  such 

points.   (For  a  graph  of  a  sample  function,  ti(x),  see  figure  2.  ) 

We  consider  a  ray  emanating  from  a  point  (0,x  )  for  which  the  group  speed, 

given  by   ^  ' — -  ,    is  negative  -  i.e.  x  decreases  as  t  increases.   From 

(15)  and  (6),  the  equation  of  such  a  ray  is  given  by 

X 

o  ^ 

(36)  t(x)     =        {  ,    .,     °  dTl. 


/^^-b2(-n) 


*     It   does  not  matter,    for  the  discussion  below,   whether       lo     is  the  parameter 
(with  X     fixed)   or   a  is  the  parameter   (with  x    =  a,    w^   =  h2(s'(a),a). 


o     o 
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At  the  turning  point,  t  =  t(x  ).   From  (36)  we  see  that  t(x  )  is  independent 
of  X  and  is  a  function  only  of  the  parameter  (w  or  cr).  Also 


(37)    ^ 
dt 


l^^Afj 


h 


For  t  >  t(x   ),    (9)  must   still  be   satisfied  and  therefore 


.A. 

(38)  t(x)    =t(xj      +    f        .^^     dr). 


X 

0) 


where  k(T],w)  is  again  a  solution  of  (6).   Because  t  must  be  a  real,  increasing 
function  of  x,  the  integral  on  the  right  side  of  (38)  must  be  real  and  positive. 
Therefore  k(T],co)  must  be  real.   This  can  only  be  true  for  x  >  x  .   Then  for 


w 


CJ 


the  integral  to  be  positive     —7 y     must  be  positive 


k(T),w) 


h 

0)  o 

Hence 


~^^^    ~     /c.^-b2(n) 


With  this   information  and   (37)^    equation   (38)   can  be  rewritten  as 


X 

"      ^  h 


(39)  t(x)      =  +  1  , ,  "  tiTi;     t>t(x). 

dx 
By  differentiating  with  respect  to  x  in  (39),  we  see  that  -r—  >  0  for  t  >  t(x  ). 

Those  rays  for  which  there  exists  a  point,  x  =  x^,  are  described  by  (38) 

for  t  <  t(x  )  and  by  (39)  for  t  >  t(x^).   The  parts  of  the  rays  described  by 


(36)  and  (39)  will  be  called  primary  rays  and  turned  rays,  respectively.   For 

^dx 
^dt 


dx 
rays  emanating  from  (0,x)  and  moving  to  the  right  (-rr  >   0),  there  is  no  turning 
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point.   In  this  case,  the  entire  ray  is  called  a  primary  ray.   These  rays 
are  shown  in  figiire  (l)  at  the  end  of  this  section.   The  contributions  to 
the  total  solution  associated  with  primary  and  turned  rays  will  be  denoted 
by  Up  and  lu  respectively. 

By  an  analysis  similar  to  that  which  was  used  to  derive  (39);  "the  phase 

of  u„  can  be  determined.   The  result  is 

T 

r  X 

O      X 


(IfO)     s-(t,x)  =  s^  :^   /  /  +  /  /^^-b^(Ti)  dTi  -  h^t  y  . 


X       X 
WO)  ^ 


In  the  case  of  one  dimensional  or  zero  dimensional  initial  data  for 
s(t,x),  the  jacobians  are  given  by  (2l)  and  (27)  respectively.   When  there 
is  a  turning  point  those  expressions  are  valid  only  for  t  <  t(x  ).   Both 
contain  an  integral  which  diverges  at  x  =  x  ,  but  in  each  case  k(x,co) 
(which  is  zero  at  x  =  x  )  is  a  multiplicative  factor  of  the  integral.   The 
limit  of  j(t(x))  as  x  -^  x  exists  and  can  be  found  by  L'Hospital's  rule. 

To  find  the  jacobian,  j(t(x))^  for  t  >  t(x  )^  (39)  must  be  implicitly 
differentiated  with  respect  to  the  parameter  (w  or  cr).   However,  if  one 
simply  differentiates  under  the  integral  sign  in  (39),  the  resulting  integral 
diverges.   Instead,  it  is  necessary  to  integrate  by  parts  with  respect  to  \\ 
in  (39),  before  differentiating.   When  this  procedure  is  carried  out  for  the 
case  of  zero  dimensional  initial  data  (with  was  parameter),  and  one  solves 
for  j^(t(x))  =  1-^  h  the  result  is 
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(l|l)      J^itix)) 


2     2,,  /2     2,    ,    (2oj  -b    (x    )) 

2u)  -b    fx)      ^    ./gj  -b   fx)    ^  ^    o' 


wb(x)b' (x) 


Jui^-\)    fx! 

CO 


cov^2-b^(x^)     b(x^)b'(x    ) 


X 


O  X 


w 


?      ? 

200  -bfn 
w,ycjj  -b    (Tj 


b(Tl)l3'U\ 


dn 


Then,  by  the  method  of  section  (1.2.2),  the  amplitude  is 


(1|2)   z(t(x))  =z 


r  J. 


-1 1/2 


JT(t(^)) 


with  jm(t(x))  and  j  given  by  (Ul )  and  (26)  respectively.   The  field  associated 
with  the  turned  rays  will  be  denoted  by  u,^.   (39)  j  (^-O)  and  f i|-2 )  describe  the 
rays,  phase  and  amplitude,  respectively,  for  u  in  the  case  of  a  zero 
dimensional  initial  manifold. 

To  find  u,^  for  the  Riemann  function,  we  use  (l+o)-(U2)  in  (32)  with  u 


+ 


+ 


replaced  by  u,^(s"  and  z"  are  given  by  (3^).)'  '^^^   result  is 

x  ^ 


cos 


(^3)   u^-v. 


U)    w 


r2nXh/co2.b2(x^)   j^(t(x)) 


1/2 


*  We  note  that  j  (t,x),  given  by  (41)  might  be  zero  for  some  (x,w). 
Equivalently,  there  are  caustics  associated  with  one  or  more  turned  rays 
in  (t,x)  space.   In  this  case,  (U3)  gives  u  on  the  turned  rays  for  times, 
t,  less  than  the  time  at  the  caustic.   By  appl^ang  the  discussion  of  phase 
change  at  a  caustic  in  [2"]  and  the  results  of  section  (2)  below,  we  find 
that  past  the  caustic,  the  phase  term,It,  appearing  in  (43)  must  be  replaced 
by  -  f  .  ^ 

k 
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u  (t,x)  is  then  given  parametrically  by  (39)  and  (U3). 

To  find  the  asymptotic  expansion  of  the  Riemann  function  at  any 
point,  (t,x),  one  must  first  determine  all  primary  and  turned  rays 
passing  through  that  point  from  (30)  and  (39).   Then  from  (35)  and 
(^3)  respectively  the  contributions  u  and  u,^  associated  with  each  of 
these  rays  must  be  determined.   The  sum  of  all  of  these  contributions  is 
the  asymptotic  expansion  of  the  Riemann  function  at  that  (t,x). 

(1.6)  An  Initial- Boundary  Value  Problem  -  The  Reflected  Field. 

We  introduce  a  boundary  at  x  =  x  and  the  boundary  condition 

a 

(kh)       u^(t,x^)  -  iX^u(t,x^)  =  0. 

The  total  solution,  to  be  determined  for  x  -  x^,  will  be  assumed  to  consist 

of  three  types  of  terms,  denoted  by  u  ,  u  and  u^.  u,^  is  that  part  of  the 

solution  associated  with  the  turned  rays  and  therefore  non-zero  only  for 

X  >  X  -  x.^.   Up  denotes  the  part  of  u  associated  with  the  primary  rays, 

some  of  which  reach  the  boundary,  u  will  be  associated  with  a  new  set  of 

R 

rays,  to  be  called  the  reflected  rays.   For  each  primary  ray  "incident" 
on  the  boundary,  a  reflected  ray  will  be  "emitted. "  It  is  assumed  that 
Up  satisfies  the  initial  conditions  and  we  therefore  require  that  the 
reflected  rays  do  not  intersect  the  initial  manifold;   if  they  did,  the 
total  field  at  (0,x)  would  be  given  by  u  +  i-l  and  would  not  satisfy  the 
initial  conditions. 

For  the  purposes  of  this  section,  it  suffices  to  assume  that 


(U5  )   u  =  Up  +  u^. 
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In  order  to  deternilne  vl^  by  the  ray  msthod,    it   is   necessary  to   find 

the  values  of  co,   k,    s  and  z   for  u     at    (t,x   ).      The  value   of  t   is   given 

by   (15)   or    (30)  with  x     =  x.^.      The  value  of  b(x)   at  the  boundary  will  be 

denoted  by  b   . 
o 

It  will  be  assumed  that  u  and  u^  are  both  of  the  forms  given  by  (2), 
The  subscripts  P  and  R  will  be  appended  on  the  phase  and  amplitude  of  u 
and  u^  respectively.   Then  from  (h^),   to  highest  order  in  \,    (1|1|)  becomes 


(1^6)   Zp(t,x^)exp[iXSp(t,Xg)](-|^  Sp(t,x^)-g  + 


Zj^(t,Xg)exp[iXSj^(t,x^)]  [-^  s^(t,x^)-^]  =  0. 


Since  (h6)   must  hold  for  all  sufficiently  large  \,    it  follows  that 


(U7)    Sp(t,Xg)  =  S^(t,Xg). 


Differentiation  of  (h7)   with  respect  to  t  yields 


(^8)      ^  ^p^^^^b)  =  -5t  ^r(^^^>^   "r  =  "p  =  -  ^OP- 


It  then  follows  from  the  dispersion  relation,  (6),   that 


(1^9)     ^(^^^r)  =  -^-p(^,^-p)- 
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If  the  signs  were  the  same  in  (^9)j  the  reflected  rays,  deternilned  as 
solutions  of  (9),  would  be  the  same  as  the  primary  rays  and  would  therefore 
intersect  the  line  (0,x).   It  has  already  been  shown  why  this  possibility 
must  be  rejected  and  therefore  the  lower  sign  must  be  chosen  in  (^9).   Then 
by  using  (5)  it  follows  that 


^^°-'    S Sc 


In  addition,  by  differentiating  (31)  with  respect  tox  (x  =  x  )  and  evaluating 


at  X  =  x^,  we  see  that 

(51)      k:(x     0))  =   ;/A^ 


and  therefore 


(52)    k-(x  c.)  =tJJ~^   ;  k-(x,a3)  =  tfKJT^) 


R'  B'  /   -  r    -^o      '        -R 


Ss 

By  substituting  the  proper  value  of  k  for  -^  in  {h6)   and  using  (^7),  it 


follows  that 


+n2.  2" 


(53)    z-(t,x^)  =  R_^z;(t,  x^)  ;   R_^  =    —=^ 


-  ^       O     ^ 


*  2   2 

The  special  case  w  =  b  requires  more  careful  investigation.   For  this 

case  k  =  k  =  0  and  the  group  speed  (^)  is  zero;  the  ray  is  tangent  to  the 

space  time  boundary.   This  case  will  not  be  considered  here  and  therefore  we 

assume  that  w  >  b^  . 
o 
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The  reflected  rays  must  again  satisfy  (9)  with  w  and  k  given  by  (ij-8) 
and  (52).   The  time  at  which  a  reflected  ray  is  "emitted"  is  the  same  as 
the  time  at  which  the  corresponding  primary  ray  is  incident  on  the  boundary. 
We  therefore  find  that 


h 


o   X 


d-T], 


where  x     is  fixed  when  w  is  the  parameter   or  x„=  o,    w  =  -   h    (s' (cr)^a), 

when   o  is  the  parameter.      With  the   same  consideration,   by  using   (lO) 

with  k     and  w     determined  above,  we  find  that 
K  K 


(55)         s^(t,x)  =  s^ 


/Ab^     dTi  -   ht\      . 


^B       ^ 


When  o   is  the  parameter  (i.e.  the  initial  data  for  s  is  one  dimensional), 
from  (22): 

1/2 


(56)   Zp(t,X^)   =   Z^[j(t(Xg))]     . 


Then  from  (53)  and  (19)   with  t-j_=  ^(^) 


(57)   ^~^=\--o 


j(t(x3)) 


-1/2  -  yt(x^-i 
yt(x)) 


1/2 
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with  jp(t(x))  calciilated  from  (5^).   It  is  easy  to  check  that 
the  jacobian  for  incident  and  reflected  rays  are  the  same  at 
X  =  x^,  so  that 


(58) 


z-  (t(x))  =R^z; 


1 

-i"2 


JR(t(x)) 


For  the  zero  dimensional  manifold  we  obtain,  in  a  similar  fashion, 
+  + 

(59)    z-(t,x)  =  R^r 


L  jj^(t(x)) 


Again,  j  (t(x))  must  be  calculated  from  (5^). 

In  the  case  of  the  Riemann  function  for  the  bounded  domain, 
X  ^  x^   ,     we  find  that 


O  X 


R     exp    {  +  ij 


(60  )     Up.--   L 


tl^- 


b       dT]   -  h  t 
o 


B         B 


+  rjr 
i"   k 


SnXhJu,  -b    (x^)      jp(t(x; 


nl 
2 


where 


(61) 


jg(t(x; 


:;2  ^2, 


X 


.  A-v"(,)   r  f   b!aii 


4 '4    C^'-A,f^' 


dri 
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as  calculated  from  ('^h).      ■LL(t,x)  is  then  given  parametrically 

ty  (5^)  and  (59).   From  (6l)  one  can  see  that  there  are  no  caustics 

associated  with  the  reflected  field  for  the  Riemann  function. 

To  find  the  asymptotic  expansion  of  the  Riemann  function 
for  the  bounded  region,  x  >  x^,  we  must  append  the  discussion  at 
the  end  of  section  (l.5)t   in  addition  to  the  primary  and  turned 
rays  passing  through  a  point  (t,x),  one  must  determine  the  reflect 
rays  and  add  to  the  solution  the  contributions,  u^,  associated  with 
these  rays.   (Reflected  rays  are  shown  in  figure  (1).  ) 

For  any  other  problem,  the  expansion  by  the  ray  method  is 
interpreted  in  exactly  the  same  way:  the  solution  for  any  (t,x) 
is  given  by  a  sum  of  contributions  associated  with  all  the  rays 
passing  through  the  point. 

The  asymptotic  expansion  of  the  Riemann  function  will  be 
obtained  in  the  next  section  by  asymptotically  expanding  the  exact 
solution.   It  will  be  seen  that  the  two  solutions  agree.   In 
subsequent  sections,  other  types  of  initial  data  and  sources  will  be 
considered,  and,  in  each  case,  the  solutions  will  agree. 

2.   Introduction  to  the  Asymptotic  Exriansion  of  the  Exact  Solution 
The  equation  to  be  studied  in  this  section  is 

(1)  u^-  u^^-  xV(x)u  =  f(t,x;X) 
with  initial  conditions 
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(2)  u(0,x)  =  u^(x;X),  u^(0,x)  =  u^(x;X) 

and  boundary  condition  given  by  (l.Uij-).   It  is  again  assumed  that  b(x) 
is  monotonically  decreasing,  b'(x)<  0;  ^(xt,)  =  b  and  lim^b(x)  =  b  . 

To  solve  (l),  we  separate  variables  by  taking  the  Fourier  transform 
of  u(t,x)  with  respect  to  t.   This  transform  will  be  called  v(x,co).   It 
will  be  shown  below  that  this  function  must  satisfy  an  ordinary  differential 
equation  in  x.   When  v(x,co)  is  known,  its  inverse  Fourier  transform 
will  yield  the  solution,  u(t,x). 

In  general,  only  the  asymptotic  expansion  of  v(x,to)  for  large 
X  can  be  determined  and  not  the  function  itself.   The  inverse  transform 
of  this  expansion  is  an  integral  which  can  be  expanded  by  standard  methods 
of  asymptotic  expansions  of  integrals  with  a  large  parameter. 

It  is  assumed  that  u(t,x)  =  0  for  t  <  0.   We  also  assiane  that 

-ct 
u(t,x)e    is  absolutely  integrable  on  0  S  t  <  oo,  for  all  c  >  0  and  that 

the  absolute  integral  is  uniformly  bounded,  with  respect  to  x.   The 

Fourier  transform,  v(x,cj),  is  given  by 

00 

(3)  v(x,w)   =      /     u(t,x)   exp[iXcot]dt,      Im  w  >  0. 

In  this  case,  the  Fourier  Inversion  Theorem  can  be  applied  to  yield 

(h)     u(t,x)  =  T~  /  v(x,cj)  exp[-iXwtJdu). 

r 

The  path  of  integration,  P,  is  the  real  w  axis,  except  in  the 
neighborhood  of  the  singularities  of  v(x,oj);  the  path  is  deformed  above 
these  singularities.   Under  suitable  conditions  on  v  we  may  rewrite  (h)   as 

00 

(5)  u(t,x)=^  ^J  v{x,t  w)   explp.\uyt]di^ 

0 
with  the  understanding  that  (5)  really  has  meaning  as  a  limiting  case  of  {k) 

when  the  deformations  around  the  singularities  of  v(x,  t  w)  shrink  onto  the 


.2k- 


singtilar  points.   Of  course,  (5)  would  "be  incorrect  if  v(x,±co)  had 
non-integrable  singularities  on  the  axis.   In  this  case  one  must  refer 
to  (k)   as  a  defining  equation  for  u(t,x). 

The  ordinary  differential  equation  satisfied  by  v(x,w)  is 
obtained  by  taking  the  Fourier  transform  of  each  term  in  (l).   The 
resulting  equation  is 

(6)  v"(x,cj)  +  X^iiZ-h^ix))   v(x,w)  =  r(x,cj;X); 
where 

00 

(7)  r(x,cj;X)  =  iXcJu  (x;X)-u  (x;X)  +  /  f (t,x;X)exp[iXwtJdt 

0 

and  ('  )  denotes  differentiation  with  respect  to  x.   When  computing 
the  transform  of  u  ,  it  is  necessary  to  integrate  by  parts  twice 

XT' 

with  respect  to  t.   Each  integration  introduces  an  endpoint  contribution 
at  t  =  0.   This  accounts  for  the  appearance  of  the  initial  data  in  r(x,w;A; 
Transformation  of  the  boundary  condition,  (l.kk),   yields 


(8)  v'(x^,w)-iA^v(x^,cj)  =  0. 

The  following  "radiation  condition"  can  be  derived: 

(9)  lim  v(x,u))  =0,  Im  w  >  0. 

X-*  00 

The  assumption  that  the  integral  of  |u(t,x)|e    is  uniformly 
bounded  in  x,  is  used  to  derive  (9).   In  addition,  it  is  necessary 
to  impose  the  restriction  that  u  ,  \i,  and  f  have  compact  support  in 
X.  As  a  consequence  of  this  restriction,  r(x,(J;  A)  has  compact 
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support  in  x.   This  will  be  asstimed  throughout  the  discussion. 

2.1  The  General  Technique  of  Solution  of  the  Time  Reduced  Equation. 
We  consider,  first,  the  homogeneous  form  of  (6): 


(10)   v"(x,w)  +  A^(a)^-b^(x))v(x,(j)  =  0. 


If  V  (x,w)  and  v_(x,w)  are  linearly  independent  solutions  of  the 
homogeneous  equation  which  satisfy  the  radiation  condition  and 
the  boundary  condition  respectively,  then  the  solution  to  (6),  (7),  (8), 
(9)  is  given  by 


X 

(11)   V(X,U))  =  -  ^^/^^   )   V^(X,U3)/    V^ii,Uy)T{i,^;^)di 

^        B 


+  Y^U,^)      /  v^(^,w)r(|,w;A)d4 


The  Wronsklan,  W(v  jV^)   =  v^v^-  v^v^  ,  is  independent  of  x. 

We  define  Region  I:   co'>  b  ,  x  >  x^.   (See  figure  (2).)  The 

O  15 

functions  v  and  Vp  are  obtained  in  this  region  by  formally  substituting 

a  solution  of  the  form  a(x,w)  exp[iXs  (x,oj)]  into  (10),  [6].   For 

b  >  CO  >  b,  ,  there  exists  a  point  x  =  x  such  that  b(x  )  =  w.   This  is 

O  1  CO  CO 


ca 


2   2   2 
lied  a  turning  point  [6].   We  define  Region  II:  b^>  w  >  b-j^,  x  >  x^ 


(See  figure  (2).  )  For  this  region,  v  and  v  are  obtained  by  the  WKB 
method  [6,7,8].   It  can  be  shown  that  v(x,co)  is  exponentially  small  for 
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real  w  except  in  Regions  I  and  II.   (The  region  were  v  is  exponentially 

small  is  shaded  in  figure  2. )  For  w  =  h  ,  x  =  x^.   Since  the  solution 

o   (jO   B 

obtained  by  the  WKB  method  is  not  valid  at  the  turning  point  the 

boiuidary  condition  cannot  be  satisfied  for  w  =  b  .   This  difficulty 

is  discussed  in  section  (9). 

In  order  to  use  (5)  to  determine  u(t,x),  it  is  necessary  to 

find  v(x,±  w).   Equation  (11)  can  be  used  for  this  purpose  if  v  (x,^  w) 

and  V  (x,±'^)  are  known.   The  expansions  of  v   {x,t^)   and  v  (x,tw) 

will  be  given  for  regions  I  and  II  in  sections  (2.2)  and  (2.3), 

respectively. 

(2.2)  v^(x,lw)  and  vJx,l^^)   in  Region  I:   w  >  b  ,  x  §  x„ 
1  ^  o B 

In  Appendix     I,    Ifert  1,    it   is   shown  that   solutions  to   (lO) 
which  satisfy  the  radiation  condition  and  boundary  condition,   respectively, 
are  given  for  w  >  b     by 


(12)     v^(x,iw)-(w2-b^(x))^  exp  [±i  A  ry^j^ -b^(Ti)     dri] 


B 


(13)     V   (x,±  w)-(u)2-b^(x))^  (exv[±Kl    y^^-h^(n)   di)   +  i/] 


B 


-expr-iA/    yu)  -b   (t)     dTi-i7]   ) ;   cot  7   =      ;4o 
^B  j  ^"  -^o 
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The  Wronskian  is 

(Ik)     W(v^,V2)/v2iXe^^^. 

Since  all  restilts  are  asymptotic  expansions,  it  would  be 
tedious  to  continually  repeat  such  expressions  as  "the  asymptotic 
expansion  of  the  solution  is  ...".   Instead,  such  expressions 
are  simply  stated  as  "the  solution  is  ...".   The  sign  of  equality 
or  asymptotic  equality  will  convey  the  proper  meaning.   This 
convention  is  adhered  to  in  the  discussion  directly  above  and  in 
all  that  follows. 

The  solutions  given  in  (12)  and  (13)  coiild  each  be  multiplied 
by  a  constant  and  still  be  solutions,  since  they  are  determined 
by  one  condition  and  satisfy  a  second  order  equation.   These  constants 
would  both  be  factors  of  ¥  and  would  have  no  effect  on  (11).   Therefore 
(12)  and  (13)  are  called  the  solutions  despite  their  non-uniqueness. 

(2.3)  V  (x,±  w)  and  Vg(x,+(^)  in  Region  II:  b^>  u)  >  b(x),  x>  x^ 


In  Appendix  I,  Ifert  2,   the  solutions  which  satisfy  the  radiation 
condition  and  the  boundary  condition,  respectively,  are  determined  by 
using  the  WKB  method.   They  are 

1 

-■J-      p     X 

(15)  v^(x,^o:)^  ^"giJ"'^^   exp  iixT /^r?(^  dTi 
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2  ,  2, 


,     +UT 


(16)  v^U,!.)^!-^^^   e^ 


exp 


ix/yco^-b^cn)  dTi  +i2 


^-   X 


X 


-exp 


-iX 


a32-b2(T)   dri-i^ 


U) 


The  Wronskian,  W  =  1. 


(2.4)  The  Asymptotic  Expansion  of  the  Riemann  Function 

As  an  example,  the  Riemann  fimction  for  (l),  on  the  domain  x  §  x  , 

B 

will  be  determined.  For  this  function  f(t,x;A)  =  0,  u  (x;A)  =  0 
and  u  (x;A)  =  B(x-x  ).   From  (8),  for  this  problem. 


(17)  r(x,tJ;A)  =  -o(x-x^). 


Consequently,  equation  (11)  reduces  to 


(18)  v(x,c.)  = 


v^(x^,w)v2(x^,w) 


^(v^jV^! 


:   x^  =  max(x,x  ),  x  .  =  min(x,x  ] 


For  w  >b  ,  using  the  results  of  section  (2.2)  in  (18)  we  obtain 
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(19)   v(x,±u))^;(2iX)"^(c.2-b2(x))^  (c^^-b^Cx^)) 


O   X 


s 


R  exp 


iX   + 


/  2   2 

'w  -b  (Tl)  dTi 


+  exp 


^B  ^B 


iiA/  yoj^-b^(Ti)dTi 


< 


R^  =  -e 


_2i7  _ o ^ 


+  /2  ^2 
^     o 


^ 


The  contributions  to  u  arising  from  (19)  can  be  calciolated  by  the 
method  of  stationary  phase  [6] •   The  contributions  from  the  terms  in- 
volving the  reflections  coefficients  R  will  be  called  u^;  the  con- 
tributions from  the  other  terms  will  be  called  u^  .   By  using  the 


method  of  stationary  phase  we  find  that  vl  is  given  parametrically  with 


parameter  t*)  by  the  equations 


R     exp 


(20)        u^~    S 


uF-h^{V\)  dT]   -   wt 


I  +  in 


-  ,     /  2  ^2,    .      / 2  ^2. 
SrtA  i/co  -b    (x)    Vw  -b    (x 


""b       ""b 


L/2      ' 


(21  )      t   = 


^B     ^B 


=-       dT];   oj  >  b      . 
b^(Tl) 


■where  the  latter  equation  is  just  a  consequence  of  the  requirement  that 
the  phase  be  stationary.   Similarly  we  find  that  u^  is  given  parametri- 


cally by 
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(22)  u^- 


(1)  ^ 


cos 


A" 


'w^-b^(Tl)  dT] 


(w  -b  (Tl): 


dTl 


< 


1/2 


^ 


(23)  t  =  J 


y'^^-b2(Ti)  ° 


< 


For  the  range  b  >  w  >  b  ,  the  functions  given  in  section 
(2.3)  are  substituted  into  (l8)  to  find  v(x,^  oj).   The 


+1  exp  -iA^  ^ ^    r/)^-b^(Tl)  dTl 


+  exp 


S 


{2h)     v(x,*  ^)'^\ 


^  xh\    J^-'o^Wi   dTl 


< 


2iA(.2_^2^x))^/Nu.2_^2(xJ)^/^ 


The  second  term  of  (2i|)  is  exactly  the  same  as  the  second  term  of 
(19),  and  therefore  the  contribution  to  u(t,x)  from  this  term  is 
given  by  (22)  and  (23)  with  one  change:  the  w  interval,  for  which- 
(24)  is  a  valid  representation  of  v,  is  the  interval,  b(x  )<  to  <  b   . 
Therefore  the  condition  which  appears  in  (23),  namely,  w  >b  ,  must 
be  replaced  by  b^>  w  >  b(x^).   A  function  u  (t,x)  can  be  introduced 
to  include  (22)  and  the  contribution  to  u(t,x)  coming  from  the  second 
term  of  (2U).   This  function  is  defined  parametrically  by 
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S 


(25)  Up  =  u^l^  t 


y^2-  b2(Tl)         < 


If  the  contributions  to  u(tjx)  from  the  first  term  of  (2^1)  are 
called  iL(t,x),  then 

^     +iexp 

pO 

(26)  u^("t,x)-Z  + 


X 

*iA   /  +  /  yw^-b^(n)  dTl-wt 


CO      (jJ 


^«i  "{     ^2   ?,  xxl/U,  2  ^2,   .^l/U 


dco. 


We  again  apply  the  method  of  stationary  phase  to  the  integral. 
It  should  be  noted  that  in  order  to  find  the  second  derivative 
of  the  phase  with  respect  to  u),   it  is  necessary  to  integrate  by 
parts  with  respect  to  t]   before  differentiating  under  the  integral 
sign.   If  this  is  not  done  the  factor  (w  -b  (t ) )^' ^will  appear  in 
the  denominator  of  the  integrand  and  cause  the  integral  to  diverge 
at  the  lower  limit  of  integration,  t)  =  x  .   The  resiiting  expression 
^°^   ^(t^x)  is  given  by  the  parametric  representation 


cos 


(27)        1^ 


oj^-b^(Tl)   dTl   -   ojt^    +  J 


<j2nAyoK?(x)    JuF-^U^) 


d  cp  ^ 

doj 


(28) 


t    = 


o         x 

X  X 

C.J 


42-b2(n) 


dri;     b  >  u)  >  b(x), 


where 
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(29) 


,2 
d  9 


doj 


2aj^-b^fx) 


2  ,2 


2u)^-b^(x  ) 
o 


2  ,2, 


b(x)b'(x>yoj'^-b^(x)    b(x^)b'(x^)yc^  -b  (x 


b(Tl)b'(]l) 


2u)^-b^fTl) 


dTl 


This  result  is  valid  only  up  to  the  first  zero  of  (29),  or 
equivalently  up  to  the  first  caustic. Past  the  caustic  the  phase 
factor  ?  in  (2?)  must  be  replaced  by  -  jj;  . 

(2.5)  Comparison  of  the  Results 

The  functions,  u^,  u  and  u^,  given  by  equations  (1.3^),  (l.i;3) 
and  (1.  6o)  agree  with  the  corresponding  restilts  of  section  2,   given 
by  equations  (2.25),  (2.2?)  and  (2.20).*  The  stationary  phase  conditions 


*  One  must  note  that  w  >  0  in  the  method  of  section  (2),   As  a  parameter, 
it  plays  the  same  role  as  h  does  in  the  ray  method. 
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of  section  2  are  seen  to  be  exactly  the  equations  of  the  various  types 
of  rays  in  section  1. 

To  find  the  Riemann  function  by  the  Ray  Method,  it  was  assumed 
that  all  primary  rays  emanated  from  the  point  x  »   This  assumption  is 
now  justified  by  the  agreement  of  the  results. 
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fl)  Rapidly  Varying:  Initial  Data  -  Ray  Method 

If  equation  (l.l)  were  originally  given  without  a  large  parameter 
A,  "stretching"  the  variables  x  and  t  by  this  factor  would  lead  to  the 
given  form  of  (l.l)  and  to  initial  data  of  the  form 

(1)  u(0,x)  =  u^fA(x-x^)],  u^(0,x)  =  Au^[A(x-x^)]. 

This  will  be  called  rapidly  varying  initial  data  f 2]  .   We  assume  that 

u  (^)  and  u  (4)  have  compact  support.   Then  the  support  of  the  initial 

data  shrinks  to  the  point  x  =  x  as  A  ->  oo.   As  in  the  case  of  the 

o 

Riemann  function,  it  is  reasonable  to  assume  that  the  rays  emanate  from 
the  point  (0,x  ).   Then  the  resiolts  of  subsections  (1.2.2,  1.3.2)  apply. 
The  rays,  phase  and  amplitude  are  the  same  as  for  the  Riemann  function, 
except  that  the  constants,  s*  in  (I.3I  )  and  z"  in  (1.2?)  must  again 
be  calculated  by  the  indirect  method  for  the  initial  data  (l).   The 
expansion  of  the  exact  solution  for  this  problem  with  b  =  constant  is 
obtained  in  [2]  and  yields  the  results 

+     /     h^   \^/2       -- 

(2)  s;  =  0,   z-  =    f a  (k Je-"^  ; 


k^  =  ^  [sgn(x-x^)]yuJ  -b  (x^)  ; 

(3)   a^(kj  =   /  [uJ4)  ±i-'  u^ii)\    e'^o^d^;  h^  =  Ao^^'^-o' 
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We  substitute  (3),  (1.26)  and  (1.28)  in  (1.2?)  and  find  that 

+    h  a+  (±  rBgn(x-x  )1  yco^-b^(x  )  )e  ^ 
(k)        z-  =   °       2-^ ° 


>      p   ,       ^  1/2 


2 


<  °  ^    (^  -b^(Tl))3/2    , 


^< 


IVom  (1.32  ),  u  (t,x)  is  given  parametrically,  with  parameter  co  ^ 
by  the  equation  of  the  rays,  (I.30  ),  and 


(5)  Up-^.  Z  z-  exp  t  iA/  /  yuF--b^{T])     dTl  -  h^t  \  +^ 


^< 


When  (3)  is  compared  with  (1.3^1),  one  observes  that  the  former 
result  is  just  the  latter  multiplied  by  +  ih  a ,  [-  [sgn(x-x  )]  s/^   -l^  (x  ] 

For  rays  in  the  direction  of  decreasing  x  this  factor  reduces  to 


(6)      p+  =  +  ^\^+r+J^^-  ^^i\ 


If  p   is  introduced  as  a  factor  in  the  simmations,   Z  ,  expressing  u 

and  vu   for  the  Riemann  function,  the  resulting  equations  are  the  corresponding 

resTolts  for  the  present  problem.  To  find  u  by  this  method,  (l.hs)   must 

T 

first  be  rewritten  as  a  sum,  L  .      The  result  is 
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(7)  u^^  L  h  a_^ 
^     T   +  o  + 


+  Jw   -b  (x  ) 


exp 


X     X 

nO      p        f 

+  -.^  I  /  ,  /   /  2  ,  2 


iA    + 


'  w 


X     X 

U)     CO 


b  (Tl)  dTl-h^t  ^  ;  ia 


|8«A  hy u)^- 


b"(x^)     4p(^(-; 


1/2 


where  j  is  given  by  (l.l|l).  The  parametric  representation  of  vl  in 

(l,6o)  is  of  the  form  il  +  ul  .   Using  these  functions  the  corresponding 
result  for  the  present  problem  is 


(8)  ^^l    v,^ 


In  (9)  CO  >  b  ,  X  >  x^  and  in  (8),  co  <  b  ,  x  >  x, 

O        ij  o  ^ 


(k)     Rapidly  Varying  Initial  Data  -  Expansion  of  the  Exact  Solution 

The  problem  introduced  in  section  (3 )  will  be  treated  by  the 
method  of  section  (2).   To  do  this,  we  calculate  the  right  side  of 
(2.1)  by  using  equations  (2.7),  (2.2)  and  (3.1).   The  result  is 

(1)   r  (x,  ±  u);  A)  =  ±  iAu)[u^(A(x-x^))  +  i  u^(A(x^-  x^))]  . 


For  u  (x)  and  U-^(x)  having  compact  support  in  x,   the  support 
of  r(x,  t  co;  A)  shrinks  to  x  as  A  -»  00.   In  this  case,  for  A  sufficiently 
large  and  x  f  x  ,  only  one  integral  in  (2.11)  is  non-zero.  Therefore 
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00 

(2)  v(x,+  0))  =  -  ^^^^^  ^J      v^(x^,t  '^)v2(x^,±u))r(^,+  c.;A)d^, 

-'■'--  00 

where 

(3)  Xg  =  X,   x^  =  4  for  X  <  x^ 

X  =  X,  Xp  =  4   for  X  >  X  . 

Using  (l),  (2)  and  (3),  v(x,;!;  w)  can  be  determined  for  the  regions  I  and 
II  defined  in  section  (2). 


(k.l)     The  Expansion  of  v(x,t^)   in  Region  I:   w>b,xi'x 

We  use  (2ol2),  (2.13),  (2clU),  (2)  and  (3)  to  write  v  as 
a  s-um  of  two  terms,  v  =  v^   +  v  ,  where 


o   X 


(U)     v^(x^-  ^)  -^   - 


R_^     exp 

t  iA  r+ 

r/o^-h^(Tl)     dn 
^B 

2iA(u)2-b2(x))^/^ 

dTl 

,         r(e. 

,-  aj;A)exp 

+iA /ycu2-b2(Ti) 

0 

2   2    l/^ 
>-  b^(4)) 


d4  ; 


X, 


(5) 


V 


(1)^  + 


rii,-   (j;A)expf+  iM  /^^-^^(Tl)  dT] 


2iA(uo^-h^(x))^/^ 


,  2  ,2..v^l/I+ 
(oj  -b  (4)) 


d|. 
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In  Appendix  II,  it  is  shown  rigorously  that  the  integral  with 
respect  to  i,,   appearing  in  (k) ,   can  he  expanded  asymptotically.  Heuristically, 

we  can  arrive  at  the  same  resiilt  by  using  the  fact  that  the  support 

2  2    ^/h  2     2     l/U 

of  r  shrinks  to  x  in  order  to  replace  (w  -b  (^))    by  (w  -b  (x  ) ) 


//  2  ?  /  2  2 

yw  -b  (r|)  dT\     by  V^  -b  (x  )  (^-x  ).   The  resulting  integral  has 


the  advantage  of  being  simply  a  Fourier  transform  of  the  initial  data. 


V  can  then  be  rewritten  as 
R 


^+'^± 


(6)    v^(x,  t   Ui)   ^ 


X    X 

o 


l/F^^^)   exp[!iA  r+  r  /K^W     dn] 


2a(uj  -b  (x)  )  '  (uj  -b  (x  )  ) 


where  a  (k  )  is  given  by  (3. 3). 

Equation  (5)  can  be  simplified  in  a  similar  manner  if  one  first 
notes  that 

(7)  x-^-Xg  =  x^-  x^  -  [sgn(x-x^)](^-x^);  x^  =  max(x,x^),  x^  =  min(x,x^). 


The  resiilting  expression  for  vl^     is 


coa 


(8)     v^^^(x,-   0))  - 


^1-    [sgn(x-x^)]ya)^-b^(x^)    J  exp[+iA    r    /co^-b^CTl)  dTl] 


2A(w  -b    (x))   '    (w  -b    (x^j) 
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(Uo2)  The  Expansion  of  v(x,j;  w)  in  Region  II:  b  >  w  >'b(x  ),  x  >  x  . 

We  use  (2.15)  and  (2.l6)  to  find  v  in  Region  II.   The  results  are 

similar  to  those  above  except  that  x^  is  replaced  by  x  and  R,   is 

(2) 
replaced  by  +i.   If  v  is  written  as  a  svia,   v  =  Vp  +  \>   then 

X 
o   X 


+i  u«_^[;yaJ^-b^(x^)      exp[+  i^   r+   r^^-^^Cn)  <iTl] 


(11)     v^(x,±c.) ^^— ^ 

^.,   2     2,    ,.l/k   I   2     2,       ..\l\ 
2A(cj  -b    (x))    '       (w  -b    (x    ))   ' 


(12)     v^2)(^^+^)   -Vp^l)(x,+  0.), 


P      ■    ■  P 


(4.^)  The  Solution 

To  find  u  ,  u  ,  and  il,  the  inverse  transforms  of  (7),  (10),  (11) 
and  (12)  must  be  expanded  asymptotically.   This  process  is  simplified 
by  using  the  results  already  obtained  for  the  Riemann  function.   If 
the  first  term  of  (2.19)  is  compared  with  (7),  we  find  that  the  latter 
co-uld  be  obtained  by  multiplying  the  former  by  p^,  where 

p  =  +  ioja  +-y^oj  -b  (x  )   .   This  factor  is  independent  of  A.  When  using 
the  method  of  stationary  phase  to  expand  the  inverse  transform  of  v  (x,*  w), 
the  factor  is  simply  evaluated  at  the  stationary  point.   Therefore,  if 
this  factor  is  introduced  into  the  sum,  2,  in  (2- 20),  the  result  will 
be  tL  for  this  problem.   This  is  exactly  the  factor  given  by  (3.7),  and  oin- 
present  multiplication  technique  is  the  same  as  was  used  in  section  (3)  to 
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obtain  vu   from  the  ray  method  results  for  the  Rieraann  function.   Since 
the  resiilts  for  u  in  the  problem  for  the  Riemann  fimction  agree, 
the  same  must  be  true  in  this  problem.   This  comparison  technique 
also  allows  us  to  conclude  that  the  results  for  iL,  and  u  must  agree 
with  the  results  of  section  (3)o 

(•3)   Oscillatory  Initial  Data  -  Ray  Method 

We  consider  now  an  initial-boundary  value  problem  for  (l.l) 
where  the  initial  data  is  designed  to  fit  the  form  of  (1.2);  that 
is 


(1)  u(x,0)/^  z^(x)exp[iAs^(x)],  u^(x,0)-' -iAh^f  s^(x),  x)z^(x)expfiAs^(x)]  ; 


h  (k,x)  =  yk^+  b^ 


(x: 


The  functions,  s  and  z  ,  are  given:  z  is  assumed  to  have  compact 
'  o      o  o 

support  bounded  away  from  x=x  .   Data  of  this  form  will  be  called 
oscillatory  initial  data  [  9] .   The  boundary  condition  is  again  given 
by  (l.i+Ji).   (For  problems  where  u(x,0)'wz  exp[iAs  "]  ,  u^(x,0)'-'Z^exp[iXs J 
with  z  ,  z  ,  s  and  s  arbitrary,  u  can  be  found  by  a  superposition  of 
solutions  of  problems  with  data  of  the  form  of  (1).   The  method  to  be 
used  in  this  more  general  case  is  described  in  [2] o  ) 

For  a  solution  of  the  form  of  (1.2),  equation  (1)  yields  the 
values  of  s  and  z  at  t  =  0: 
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(2)  s(0,a)  =  s^(a);  z(0,a)  =  2^(0). 


This  eqaiation  gives  s(0,x)  initially  on  the  one  dimensional  manifold, 
t  =  0,  so  that  we  may  use  the  results  of  sections  (l.2ol)  and  (I.3.I) 
to  obtain  iL(t,x).   The  resulting  parametric  representation  (with 
parameter  a )   is 


(3)  u^^  z^(a)[j(t(x))]"^/2  exp 


iA  ^  s  (a)   + 


\; 


S 


< 


^cJ^-h^(Tl)  ^"H  -  h^t 


S 


(h)     t  = 


/2  ,2, 


dTl,  X  =  min(x,a),  x.  =  max(x,CT). 


x^  y^  -h  (Ti: 


j(t(x))  is  given  by  (1.21 ).  To  insure  that  there  are  no  caustics  we 

assume  that  s"  >  0.  The  solution  will  contain  the  features  of 
o 

reflected  and  tvirned  rays  if  it  is  also  assumed  that  s'(a)  <  0.   In 
this  case,  from  (1.9),  (1.13)  and  (l.lij-) 


(5) 


dx 
dt 


[sgn  s^(a)]  yu)  -b  (a] 


<  0;  w  =  h  (s' ia),o) 


=   ys'2(a)+b^(a: 


(0,a) 
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2   2 
For  cj  <  b  ,  each  primary  ray  gives  rise  to  a  turned  ray.   From 

(1.39  ),  the  equation  of  the  tiArned  rays  is 


(6)  t(x) 


O  X 


dT\;   t  >  t(x  ). 


X   X   yco^-b^(ri) 


The  phase  on  the  turned  rays  is  given  by  (1.^0),  where  the  choice 
of  sign,  (t),  refers  to  sgn  w.   For  the  present  case  sgn  oj  =  1. 
The  amplitude  on  the  turned  rays  is  given  by  (1.22  )  with  the  jacobian, 
j  (t(x))  calciilated  from  (6)  by  the  method  described  in  section  (1.5). 
The  result  is 


(7)   j^(t(x)) 


yuj  -b  (x) 


/  2  -^2,  s 
'w  -b  [a) 


I 


0)  -b  (x) 


(s^(a)s;;(a)+b(a)b'(a)) 


(jj 


2co^-b^(a) 


2w^-b^(x) 


b(a)b'(a)yco^-b^(a)      b(x)b' (x)yw^-b^(x; 


V. 


cr      X 

X        X 
(jj        to 


\   2u)^-b^(Tl) 


-> 


1^(11)^(11)/  y^^--b2(Tl) 


dTlS 
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With  these  results,  u  (t,x)  is  given  parametrically  by  (6)  and 


(8)  u^^zJa)[j^(t(x))]-l/2  exp 


iA  (  sja)  + 


0              X 

^    - 

-  h  t  > 
0  / 

J   ^ J  ^^"^^^^^  ^'^  ■ 

X             X 

CJ               CO 

^  -^ 

The  reflected  rays  are  given  hy  (1.5^)  with  x  =  a: 


a  X 

(9)  t  =  [  ^     [ 


h 


dT], 


Xg    x^  Jui   -b  (Ti; 


The  phase  for  tjl  is  obtained  from  (I.55)  with  sgn  w  =  1.   z(t(x))  is 

given  by  (I.58).  The  jacobian,  j  (t(x))  =  l"^!,  "i^^t  be  calciaated 
from  (9).   The  result  is 


(10)   jp(t(x)) 


^J^-b^fx) 


y. 


co2-b2(a) 


■^^^^\^[s;(a)s^(a)  +b(a)b«(a: 


a  X 


b^fTi: 


x^    x^  (oo2.b2(Ti))3/2 


dTl 


These  restilts  allow  us  to  express  Uj3(t(x))  parametrically  by  (9)  and 


-kh~ 


(11)  u^-  R_^zja)[jp(t(x))]"^/2^xp 


a       X 


(o)  +  /  +  /  JuF-h^Cn)     dTl  -  h  t 


^B  ^B 


with  R  given  by  (I.53  ). 


To  find  u(t,x),  we  add  all  contributions,  u  ,  il  and  ii  ,  on  rays  passing 
throiigh  the  given  point  (t,x).   The  parameter  w  in  (3),  (10)  or  (ik) 
is  determined  ftom  the  equations  of  the  various  types  of  rays,  (h),    (7) 
or  (10). 

(6)      Oscillatory  Initial  Data  -  Expansion  of  the  Exact  Solution 

(6.1)  The  Method  of  Solution 

In  this  section  we  shall  show  how  the  problem  introduced  in 
section  (5)  is  solved  by  the  method  of  section  (2).  The  function 
v(x,  -  to)  is  given  by  (2.11)  with  r(x,i'  w;  A)  determined  by  (2.6) 
and  (5.1): 


(1)   r(x,  +  u);  A) 


,^L.^^i.'^M..^.^M.M.^s^M]. 


The  only  new  feature  in  our  method  of  solution  will  be  that 
u  ,  iL  and  iL  must  now  be  determined  by  the  method  of  two  dimensional 
stationary  phase  [2] .   The  details  of  calculation  will  be  carried 
out  only  for  the  function  u.^.  u  and  u,^  can  then  be  computed  in  a 
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similar  manner.  As  in  o-ur  previous  examples,  special  care  must  be 
taken  in  finding  the  second  derivative  of  the  phase  function  associated 
with  iL.   This  problem  was  discussed  for  the  Rieraann  function  and 
is  overcome  in  the  same  manner  as  was  used  in  that  case. 


(6.2)  The  Function  v. 


R 


To 


find  u^  we  first  look  for  that  part  of  v  which  depends  on 


R,  .   This  function  will  be  called  v"  and  it  appears  as  part  of  the 
i  R 

solution,  v,  for  Region  I:  w  >  b  ,  x  ^  x  .   It  is  given  by  the 


equation 


(2)   v^(x,  t  cu)^  + 


x^ 


R  exp 


co^-b^(Tl)  dT]  >  +  iAs^(l) 


2(w  -b  (x))  '   (u)  -b  (I)) 


.  (^^  u)  +  hjs;(^),  4)j  z^(e)exp[iAsJ^)]de  , 


with  R    given  by  (I.53). 


(6.3)  u^(t,x; 


With  the  above  solution  for  v  ,  we  use  (2.1|)  to  define 


A 


3)     u^(t,x)^    +     2n      /       ^R*-^^    -   '^^   exp[+iAwt]dcJ. 


-1+6- 


The  phase  functions  of  the  integrands  in  (3)  are 


(h)       cp^  (o),  I) 


.fr. 


^oj^-b^(ri)  dTl-wt 


^sj^)  , 


^B   ^B 


with  first  derivatives  given  hy 


(5)    ^     =     '- 


+ 


^ 


^ 


B    B  y^j^ -13^(11) 


dTl  -  t  ); 


Scp+ 


2  ,2, 


(6)  ^  =   +yu.  -b'^(4)   +  s^(|). 


Setting  these  derivatives  equal  to  zero  yields  the  equations 


X 


(7)   t 


x„  "B  yoJ^-b^(Tl) 


dll;    s^(e)  =  +ya)^-b^(e) 


B 


Since  we  have  assumed  that  s'(^)  <  0,  there  can  only  be  a  stationary 

o 

point  for  the  upper  sign  in  (7).   Also,  by  solving  for  w  we  obtain 


(8)  u)  =  ys'2(o   +b2(|: 


=  hjs;(|),|). 
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In  addition  to  (7)  and  (8),  it  is  necessary  to  find  the  determinant 


of  the  matrix  of  second  derivatives  of   cp_|_. 


>2    x2      f-^2      > 


:^t2   ^  2 


.acua^y 


The  matrix  will  be  called  S  9+  .   At  the  stationary  point  defined  by 
(7)  and  (8),  we  find,  by  using  (5)  and  {G) ,   that 


(9)  det(a  cp_^) 


2  •k2/£' 
w  -b  (^, 


b(^)b'(g)+  s^^{i)sl{i)f 


M-^^{i) 


2  ,  2,__  ,,3/2   '" 


B   B 


(w  -b  (n))- 


The  first  term  of  (9)  is  clearly  negative  for  co  >  b  .   As  in  section  (5) 

we  assume  that  s'(4)  <  0,  s"(4)  >  0  and  throughout  our  discussion 

2 
b(^)  >  0,  b'(|)  <  0.   Therefore  det(a  ^^)  <   0.   This  implies  that 

2  2 

the  eigenvalues  of  (S  cp  )  are  of  opposite  sign  and  hence  sig(a  '^   )   =0. 

Using  our  results,  (2)  (9)  and  the  formula  for  n-dimensional 

stationary  phase  in  [2]  we  find  that  iL(t,x)  is  given  parametrically 

by  (7),  (8)  and 


X 


h^(s;(^),|)R^z^(|)exp 


iA  sj^) 


+   /  +  /  yaj^-b^(Tl)  dTJ-urt 


10 


V 


B   B 


2  ^,2,  , 
w  -b  (x) 


^.^-b^(^: 


det(a  cp^) 


1/2 
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{6.h)      Comparison  of  Results 

Comparison  of  (7) (10)  vd.th  the  solution  for  u^(t,x)  obtained 

by  the  ray  method  reveals  that  they  agree.  The  calculations  of  Up 
and  TJL  have  also  been  carried  out  and  we  find  that  they  also  agree 
with  the  resiilts  obtained  by  the  ray  method. 

(7)  The  Green's  Function  -  Ray  Method 

We  consider  now,  the  inhomogeneous  equation 

(1)  u  -  u,  ,-  aV(x)u  =   6(x-x^)5(t-t  ) 
^       XX   tt  o      o 

with  initial  data 

(2)  u(t^,x)  =  u^(t^,x)  =0;  t^  <  t^. 

At  X  =  x^  the  boundary  condition  (l.i|i|  )  is  imposed. 
B 

We  ass-ume  that  the  primary  rays  emanate  from  the  point  (t  ,  x  ) 
and  for  each  primary  ray  incident  on  the  boundary  a  reflected  ray  is 
produced.   Also,  tirrned  rays  will  arise  as  described  in  section  (1.5^ 
If  we  set 

+ 

(3)  Up  =  u_^  +  u_;  u_^^  z"  exp[iAs-]  , 


then  the  fwictions  s"  and  z"  are  determined  by  the  methods  of  sections 
(1.2.2)  and  (1.3.2)  except  for  the  constants  s-  and  z  .   These  numbers 


can  be  determined  by  the  indirect  method. 
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(1.30)  is  the  equation  of  the  primary  rays  emanating  froni  the  point 
(0,x^).   If  t  is  replaced  by  (t-t  ),   we  obtain  the  equation  of  the  primary 
rays  emanating  from  the  point,  (t  ,x  ).   In  fact,  if  this  substitution 
is  made  throughout  the  development  of  the  solution  for  the  Riemann 
function,  the  corresponding  components  of  the  Green's  function  result. 

For  example,  from  equation  (  1.33  )^  u  (t,x)  for  the  case  of 
constant  coefficients  is  given  parametrically  by 


+      -1/2 
(1|)  Up-  Z  z-(t-t^)     exp 


iAJ  s^t^/^l^'lx-xj-hjt-tjl 


x-x  \-   ^ — (t-t  ). 

o'     h     ^   o 
o 


The  expansion  of  the  exact  solution  for  b  =  constant  has  been  carried 

out  in  [2].   Section  (h) ,    equation  (5)  of  that  paper,  with  a(k,l  h  )  =  1 

2  -3 
and  h"(k)  =  b  h  ,  yields  the  proper  solution  for  comparison  with  (h). 

+  ~+ 

s-  and  z"  are  then  seen  to  be: 
o 


+ 


.      /     h     \^/%^" 


(5)    s-    =   0,    r 


\ 


0_ 

8jrAb^(x 
^  o 


Comparing  (  5 )  with  (I.32  )  we  see  that  they  differ  only  in  the 
factors  e"*"-"-"  in  the  equations  for  z*.   It  therefore  follows  that  the 

Green's  function  can  be  extracted  from  our  resaalts  for  the  Riemann 

+ 
function  by  replacing  t  by  (t-t  )  and  multiplying  z"  in  the  restilts 

for  the  latter  by  e"*" 
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8.   The  Green's  Function  -  Expansion  of  the  Exact  Solution 

In  this  section,  we  consider  the  asymptotic  expansion  of  the  exact 
solution  for  equation  (7.1)  with 


(1)  A^f(t,x;A)  =  5(x-x  )6(t-t  ), 


homogeneous  initial  data,  and  the  "boundary  condition  defined  ty  (l.lj.i|), 
It  is  now  necessary  to  solve  equation  (2.6),  where  r{x,t   w;  A), 
determined  from  (l)  and  (2.7),  is  given  "by  the  equation 


(2)  r{x,i   u;A)  =  6(x-:^)e"^^'^o 


Ey  comparing  (2)  with  equation  (2„l8),  we  see  that  our  solutions  for 
v(x,-  w)  must  be  just  multiples  of  the  solutions  for  v  in  the  expansion 
of  the  Riemann  function.   The  multipliers  are 

-exp[-iAwt  ]  =  expf^  iit]  expFj:  iAurt  J.   If  this  factor  is  inserted  in 
the  integrals  in  equations  (2.20)  and  (2.26  ),  the  functions  thus 
obtained  will  be  the  analogous   expressions  for  the  Green's  function. 
The  constant,  exp[+  inj ,  is  just  a  multiplier  of  the  terms  in  the 
summations,  Z  ,  in  (2.20)  and  (2.26).  The  factor,  exp[-  iAurt  ],  has  the 
effect  of  replacing  t  by  t-t  in  the  integrals  used  to  determine  the 
Riemann  function.   Therefore,  the  same  adjustment  must  be  made  in  the 
stationary  phase  conditions  and  parametric  representations  of  the  solution. 
These  are  the  same  observations  as  were  made  in  section  (7)  to  calculate 
the  expansion  of  the  solution  by  the  ray  method.   It  therefore  follows 
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that  since  the  two  methods  led  to  the  same  solution  for  the  Riemann 
function,  the  same  must  be  true  for  the  Green's  function. 
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9.   Diffraction  in  Space -Time 

When  finding  the  solution  to  initial -boundary  value  problems  by 
the  methods  introduced  in  sections  (l)  and  (2),  all  rays  were  confined 
to  that  portion  of  (t,x)  space  where 


X 

X     O      - 

b 

dT\. 


/ 


Xp   A  2  ^2, 


^B  ""B  yb  -b  (ri) 
o 

This  can  be  seen  by  examining  figure  (l)  or  the  equations  of  the  various 
types  of  rays  in  sections  (l)  -  (8). 

In  addition,  we  noted  in  section  (2)  that  use  of  the  WKB  method 
to  solve  the  time-reduced  form  of  (l.l)  had  the  disadvantage  of  forcing 
us  to  exclude  the  values  co  =  ±bo  from  consideration.   Both  of  these  dif- 
ficulties can  be  overcome  simultaneously  by  a  method  used  in  [3]  to  treat 
an  analogous  problem  for  the  reduced  wave  equation.   We  note  that  the  ray 
defined  by 

X 

(2)        t  (x)  =     ^--=% dil,  t        <  t(xj; 


Ao-b'(^) 


<  -  ^  B' 


X        o 


(3) 


^>^^)=f^  \  /2%  dri,  t>t(x^). 


is  tangent  to  the  space  time  boundary.   This  ray  will  be  referred  to  as 

the  critically  reflected  ray.   Equivalently,  in  space,  this  ray  is  incident 

dx 
on  the  boundary  at  velocity  IT  =  0.   We  shall  show  how  the  Riemann  function 

is  determined  for  t  >  t  (x). 

We  consider  the  initial -boundary  value  problem  for  u  defined  by 

the  equations  (l.l),  (1-29)  and  (l.i+i+).   We  define  v(x,co),  the  Fourier 
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transform  of  u,  by  (2.3)  and  its  inverse  transform  by  {2.k).      Then,  as 
shown  in  section  (2) ,  v  is  a  solution  of  the  problem  defined  by  the 
equations  (2.6)^  (2.8)  and  (2.9)  with 


( ^4)       r(x,wjA)  =  5(x- 


■X 


o 


As  in  section  (2)   we  require  that  v  satisfy  the  radiation  conditi 
(5)       lim  v(x,(jj)  =  0,  Im  w  >  0. 

X->  00 

It  is  shown  in  section  (2)  that 

ff.^  v^(^,c.)vg(x^,u))  .   ,     , 

(,Dj       V  =  :  X.  =  max  (x.x  )  ,  x^  =  min  (x.x  )  ; 

W(v^,V2)        ^  °    <         '  °  ' 

(7)  WCv^^v^)  =  v^v^-  v^v^; 

where  v  (or  v  )  is  a  solution  of  the  homogeneous  form  of  (7)  satisfying 
the  radiation  (or  boundary)  condition.   Our  goal  is  to  obtain  asymptotic 
expansions  of  v  and  v  ( and  therefore  of  v)  which  are  uniformly  valid 
through  the  turning  point,  x  :b(x  )  =  oj.   The  method  of  obtaining  these 
expansions  is  discussed  in  [3]  and  [5]. 
We  define 

X 

(8)  I  [cp(x,aj)]3/2  =  r  (^2_  ^2^^^^i  ^^  ^  t(x,a));  cu  in  D:  b^<  oo^<  b^ 


on 


X 


and  assume  that  tir(x,u))  can  be  analytically  continued  in  the  complex 
w  plane.   For  fixed  x,  it  is  shown  in  Appendix  III  that  t   has  a  branch 

point  at  w  =  ±  b(x).   We  take  the  branch  cut  to  run  from  -  b(x)  to 
b(x)  and  define  t  to  be  positive  for  b  >  w  >  b(x) .   From  [3],  we 
define 


(9)  Y^(x,co)    =   [cp'(x,u))]"i  Ai[A^/3  e"%(x,co)], 

JTT 

(10)  Y^(x,u))    =   [cp'(x,w)]"%i[A^/3  g      3     q)(x^co)], 
Ai   is  the  Airy  function  defined  by 

(11)  Ai(a)    =     ~     \    cos  ( ~r  +  at)  dt;   a  real. 


Y  and  Y  are  asymptotic  to  two  linearly  independent  solutions  of  the 
homogeneous  equation  (2.10);  see  [6].   It  can  be  shown  that  for  Re  w   in 
D  and  0  <  Im  w  <  <  1,  Y,  satisfies  the  radiation  condition  (5).   Therefore, 
we  conclude  that 

(12)  v^(x,a))  ~  Y^(x,co). 

Let  us  introduce  the  "boundary  operator"  fi  by  the  equation 

(13)  fiv(x,ca)  =  v'(xg,u))  -  iA^  v(Xg,u)). 
It  then  follows  that 

(lU)       V2(x,co)  ~  Y^  n  Y^-Y^  n  Y^. 

Using  the  Wronskian  of  Airy  functions  given  in  [6]  we  find  that 

ijt 

(15)  W(v^,  Vg)  =  {2n)-\n   Y^)  A^/^  ^  ^  ^ 

ajid  then  (6)    becomes 

irt 

(16)  v(x,u^)~  2«A"2/3  g  Y^(x^,w)    Y^(x^,u))-Y^(x<,u)){fiYo}{fiYi} 

Using  (9)  and  (lO),  it  can  be  shown  that 

fiY         a  Ai[A^/3  e"^%(x,co)] 
(17) 


2/^ j^ , 
nY-^       fi  Ai[A  ' -^  e' -^  cp(x,w)] 
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Let  us  assume  for  the  moment  that  the  contour,  P,  in  (2.i+)  can 
be  "closed"  in  the  lower  half  of  co  plane  for  t  >  t(x  )  .   It  would  then 
follow  that 

(18)  u(t,x)  =  A(2n)"  \  v(x,w)exp[iAcot]dk'  +  2  Ai  R^^, 

where  B  is  an  integral  around  branch  cuts  and  on  a  semi-circular  path  in 

the  lower  half  of  the  w  plane.   The  functions  PL  are  the  residues  of  (18) 

at  the  zeroes  of  the  denominator  of  (17).   Using  (9),  (lO),  (16)  and  (l?), 

we   find  that 

l__ijt  i.     _  ilL'  2  !£ 

2nA3  e      ^     Ai[  A^   e     ^  cp(x,ca  )  ]Ai[  A^  e     ^     cp(x   ,m   )] 

(19)  AiR^~   -^ ^ ^-^^  exp  [-iAc^t] 

2 


From  the   definition  of  fi,   we   find  that  w     must   satisfy     the   equation 

1 

Ai(r  e"^^)  a3   te^"^/^ 

HI  ^ 


(2°)  Ai(r  e-^'^)      '     cp'(x^,.   ) 

m  ^    B'   m' 


(21)  r     =   a2/3  e2«^/3  cp(x^,(.   ). 

^  m  '    B'   m' 

It  can  be  verified  that  for  each  oj  such  that  Rew  >  0, 

m  m 

its  image,  obtained  by  reflection  through  the  imaginary  axis, 

Rew  ->  -Reco  ,  is  also  a  root.   We  will  therefore  consider  only  those 
mm'  ■^ 

roots  where  Rew  >  0  and  indicate  later  how  the  results  are  to  be 

m 

adjusted  to  include  residues  R  for  which  Reco  <  0. 

m  m 

Let  us  define  x  =  x  .   From  (8)  and  (21)  we  see  that  for  r 

m 
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to  be  finite  as  X  ^»,  x^  ~  Xg.   In  this  case  cp(xg,  w^)  and  ^(x^,  co  ) 
can  be  approximated  by  integrating  by  parts  in  (8).  By  this  process 
it  can  be  shown  that 

(22)  (.  ~  b  +  r  (2b  )-^/3(b'A-^)2/3  e"^^/3; 
^      '  m    o    m^   o^      o   '  ' 

(23)  b^  =  -b'(xp  >  0. 

It  can  also  be  shown  that  asymptotically  (20)  becomes 
Ai(r  e   )      Ce   ' 

III 

^''^       Ai(r  e-^'^)  ^  (2b  b'A-l)V3   ' 
^  m    '    ^   o  o 

In  addition,  we  may  simplify  the  denominator  of  (19)  by  noting  that 
(25)    -t  "Ai(  V"-)  ~  2V3(,^,J,-l)-2/3  ,^,^i(, 


The  Airy  functions  appearing  in  (19)  can  be  expanded  asymptotically  by 
using  the  results  in  [6].   We  then  find  that 

/2bb'\2/3e'^V6^^i(^3-5«i/3) 

( 26)   AiR  ~  -^    ^tt; 

"^  \  W      2b  fiAi(r  e  "^) 


1'-       ^     ^^  x. 

y  +  y  ((.2.  b2(Ti))*  dTi-a>^t| 


[q)^(x,u3j9'(x,c.Jcp^(x^,,.Jcp'(x^,c.J]^ 


The  denominator  in  (26)  can  be  further  simplified  by  differentiating  (8) 
and  evaluating  the  result  at  x  and  x^.   In  addition,  the  exponent  in  (26) 
can  be  expanded  asymptotically  by  the  sajne  procedure  as  appears  in  [3]- 
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XX  ,-    X        X 

o  ^  I  o 


(27)  iA<    r       +\        (CO^-    b^(Tl))2dTl-ujt)~     iA<    C+r       (b^-    b^(Tl))idTl     -    b    t 


-j-o)"/^(^;)^/^^'/^e^"/'(vv-T^/' 


where 


X  X 

(28)  t      =      C  ,  g°      ,  dr,    ;    t    -  T     .    r  ,      °  d,. 


t      is  the   quantity  previously  defined  as  t   (x   ),   the   arrival  time   of  the 

critically  reflected  ray  at  the  boundary.      We   interpret     T^  as   a  time  of 

B 

emission  of  a  new  type  of  ray  which  we  will  call  the  time -diffracted  ray. 

The  requirement  that  t  >  t  (x)  implies  that  T  >  t   and  we  see  that  the 

->  B    B 

time -diffracted  ray  is  just  the  critically  reflected  ray  displaced  in 
space-time  by  the  amount  T^-  t  .   (See  figure  (l)). 

D  B 


We  may  now  rewrite  (26)  as 


,   ^2/3     ^ni/6 
2b  b      ^ 
(29)        AiR       '         °  ° 


'fiAi(r^e"5'''-/3)exp  iAJ£+  ^/b^-  b2(  t])  drj-b^t  ( 


m     \        A 


2b  fiAi(r  e"^'^)[b^-b^(x)]2[b^-b^(x   )]2 
o  m  o        ^'o  o 


exp  r  (Ab'2[2b    VY^ ^   e^^^^^^T   -t    )-^r^/^] 
^l^m^      o   ^      o-"      ''  ^-"b     B^       3       m    J 
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For  the  special  case  where  ^  -^  «  in  (2U) ,  r  >  0  [6] .   In  this 


m 


.1/3. 


case  from  (29)    R     decays   exponentially  with  A  '  ^rJ,T^-  t^) .      A  further 
simplification  is  possible   if  we  assume  that   r^  =   0(A  ) ;    0  <  a  <   1.      In 
this   case  we   can  use  the  Wronskian  of  the  Airy  functions   [6]    and  the 
boundary  condition  (I3)    to   simplify  the   quotient   of  boundary  operators 
in  (29).      The  result   is 


(30)        "^^^V 


-5Ki/3x 


^2/3^-iV6 


fiA  i(r^e        ) 


2rt(2b  b') 
^      o  o 


1/3 


Ai(r  e""^'^)flA'i(r  e   '^^ 
m  m 


By  using  (13,    (2U)    and  the   asymptotic   form  of  (p' {x^,ujj  ,   we   simplify  this 
expression  still  further  and  obtain 


fiAi(r  e"^''^/^)  . 

(31)  2L_: |2ne'^^/2 


fiA'i(r  e 

m 


■17t, 


2b  b')2/3r  e^^/W/3^2" 


00  m 


AL2(r   e"^'^) 
m 


Substitution  of   ( 3I)    into   (29)    yields,    for  large   values   of  r   ,   the   expression 


(32)  XiR  ~ 


U  y.\    2/3 


-iJt/3 


2b  b 


exp 


m  X 


o  o 


b^-    b^(Tl)     dTl-b^t 


it-jtb 


b^-  b^(x) 
o  ^ 


iW 


b^-  b^(x   ) 
o  o 


lA 


a^(r^(xV^[2bJ-^)V3  , 3.1/6,,^.,^) 


3       m        f. 


[2bb')2/3r  e^i/3+x2/V' 
^00  m  ^  _ 


..2/    -xus 
Ai  (r  e   ) 
^  m 


We  show  in  appendix  III  that  the  branch  cut  integrals  and  the  integrals 
along  the  semi -circular  path  in  the  lover  half  plane  do  not  contribute 
to  the  result  and  therefore  from  (22) 
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(33)      u(t,x)~  Z  AiR^ 
m=l 


where  R^  is  given  by  (29)  for  r  =  0(l)  and  by  (32)  for  r  =  0(  A*^)  > 

0  <  a  <  1. 

To  find  R  for  Rew  <  0,  ve  introduce  the  following  rotation 

(3^)      ^^[f  (W3_  +  iw^)]   =  f(-Lo^  +  ia)^). 

It  ean  then  be  shown  that 

_  2ni 
(35)      V(^iRj  =^^   ^  '        %^«'- 

Application  of  (35)  to  (32)  allows  us  to  find  the  residues  for  Reoj  <  0. 


m 


For  r   =  0(1),  we  replace 


UAifr  e         ^     )  by  OAilr     e       ^    j      in  (29) 


m       /  \  m 


V' 


and  call  the  new  expression  R  .   It  can  be  shown  that 

m 


1  TT 

3 

m'      m        '     m 


(36)       */   (AIR  )  =  AIR  e   3   ;    r  =  0(l), 
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Appendix  I 

The  purpose  of  this  appendix  is  to  find  two  linearly  independent 
solutions  to  the  homogeneous  ordinary  differential  equation,  (2.10).   One 
solution,  V  (x,±co),  should  satisfy  the  boundary  condition,  (2.8),  and  the 
other,  V  (x,±(jj),  should  satisfy  the  radiation  condition,  lim  v(x,co)=  0, 

X  -*oo 

Im  CO  >  0. 


Part  1:   v  (x,±u])  ,  v  (x,±w  )  for  tj  >  b  . 

By  formal  substitution  into  ( 10)  of  a  function  in  the  form 

a(x,u))exp  iX[s(x,(jo)],  one  finds  that  two  linearly  independent  solutions 

have  asymptotic  expansions  given  by 

_  1        ^ 

(1)  v^^^x,co)  =  (co2-b2(x))   ^exp[ixC  {u>^-h^iv,)^  d^] , 

_  1  X 

(2)  v'^^^(x,u))  =  (oj^-b^(x))   ^expL-ixV  (co^- b  ^(  t])  2  dt]] . 

To  find  a  solution  satisfying  the  radiation  condition,  (2.9),  we 
must  extend  v^    and  v^    into  the  upper  half  co  plane  and  examine  them 
for  Im  w  >  0  and  x  ^  00. 

For  complex  00,  the  asymptotic  behavior  of  v    and  v  '    depend 
on  the  nature  of  the  exponents.   We  define 

X  1 


(3)      e^(u))  =^^  i^^-   b^(^)^  dri. 


^B 

2    2    — 
The  function,  (10   -  b  (ti))2^  is  defined  as  an  analytic  function  of  w  in 

the  "cut  plane",  with  branch  cut  from  -b(T])  to  b(Ti).   Also 
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r 


(h) 


(c.2.  b2(Ti))2  =y 


/  2    2 

'CO  -  b  ( ri)  ,  u)  >  b(Ti) 


/  2   2 

'co  -  b  (r|)  ,  -  uj  <  b(Ti) 


"^ 


For  cu  =  u.'  +  i  CO  and  co  >  0  but  small,  e^(aj)  can  be  approximated  by 


(5) 


e,(co)  z   e,(co  )  +  e'(co  )  (oj-oo  ) 


±~  O' 


With  this  approximation,  we  find  that 

(6)       Im  e  (co)>  0,  Im  e  (w)  <  0;  Im  co  >  0. 

We  conclude  that  v   (x,co)  decays  exponentially  as  x  -^  <»  and  v   (x,co) 

2   2  — 
grows  exponentially  as  x  ->  oo.   By  using  the  definition  of  (cj  -  b  )  2  in 

7  2   2 
CO  -  b   we  rewrite  (l)  and  (2)  as 


(7; 


V   (x,±co)  ~  (co  -  b  (x))     exp 


±iA  \   /co^-  b^(Ti)  dT] 


(8) 


(2),   ,  ,   ,  2  ?,  .  "  ^ 
V   (x,±cli)~  (co  -b  (xj)    exp 


iA  \   /co  -  b  ( rj)  dT\ 


J  CO  >  b 


The  conclusion  of  the  analysis  above  is  that  the  function, 


(9) 


;^(x,±co)  =  v^^^(x,±co), 


satisfies  the  radiation  condition  ,   In  (7)  and  (8),  one  should  also 

2    2     "5" 
define  the  factor  (co  -  b  (x))    in  terms  of  the  positive  fourth  root, 

/CO  -  b  .  Any  phase  factors  introduced  by  this  definition  will  also 


appear  in  the  Wronskian  and  would  not  affect  v(x,cj)  as  defined  by  (2.11) 
We  therefore  interpret  the  amplitudes  as  positive  fourth  roots. 

The  function,  v  (x,±co),  satisfying  the  boundary  condition,  is 
given  by  the  following  linear  combination  of  v    and  v   : 


(10)     yJx,±u^)    =   ± 


v^"'"^(x,±  co)exp[±i7]-v^^\x,±u)exp[+  ij] 


,  CO  >  b 
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7  is  determined  by  substituting  (lO)  into  the  boundary  condition,  (2.8) 
For  this  purpose,  the  following  form  of  (lO)  is  simpler  to  use: 


_  1 
(11)     Y^(x,±ijj)~   2i(aj^-  b^(x))    sin 


\  /cj  -  b  (ri)  dTi  +  7 


7  is  then  determined  by  the  equation. 


(12)     cot  7  =  /^   g   ;  b^  =  b(x^) 


/ 


w  -  b" 


The  Wronskian,  v  v '  -  v  v'   can  be  obtained  by  using  (9)  and  (lO) 
The  result  is 

(13)     W(v^,V2)~  2iAexp[;i7]. 


Part  2:   v^(x,±(jo),  v^(x  +co)  for  b>co>bT 
1^  '    '   2   '         o     1 

For  each  co  in  this  range,  there  exists  a  turning  point,  that  is 

2   2  2   2 

an  X  =  X  ,  such  that  lo  -  b  (x)  >  0  for  x  >  x  and  uj  -  b  (x)  <  0  for 

co'  CO 

x  <  X   .      V   (x,±co  )    is  determined  by  its  behavior  for  x  large.      In  partic- 

U)     1 

ular,  for  each  fixed  co,  we  can  consider  the  range  x  >  x  to  determine 

'  '  CO 

this  function.   The  analysis  of  part  (l)  is  simply  repeated  with  x 
replaced  by  x  and  the  result  for  x  >  x  (analogous  to  (7)  and  (9)  )  is 

■^  CO  CO 

_  1 

2   2    ~   h 
(lij-)     V  (x,±(io)  ~  (co  -  b  (x))     exp 


:iA  \  y cj  -  b  ( Ti)  dT)  ; 

CO 


X  >  x 


V  (x,±co)  is  determined  by  its  behavior  at  x  =  x^  <  x  ,  where 


cj  -  b  (x)  <  0.   Analogous  to  equation  (ll)  we  may  write 

_  1  X 

o       o  ~ 

(15) 


2       2        h 
V  (x,  CO  )~  2(b  (x)  -  CO  )    sinh 


A 


I,/^ 


ll)  -  CO   +7 


B 


By  using  the  boundary  condition,  (2.8),  we  find  that 
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(16) 


coth  7  = 


/ 


"b  -  CO 


In  order  to  determine  v  (x,±  w)    for  x  >  x  we  must  use  the  WKB  connection 
formulas  [6,8].   For  a  homogeneous  second  order  ordinary  differential 
equation,  these  formulas  are  asymptotic  expansions  of  a  specific  pair  of 
linearly  independent  solutions  for  x  <  x  and  x  >  x  .   The  connection 
formulas  for  (2.i0)  are  given  below. 

1        X  , 


(17) 


v^3)(^)^(^2_  ^2^^^^  k   ^.^ 


] 


A  \  /co2_  ^2^^^  ^^  ^  I 


'  ""^   \' 


(18) 


v^3)(^)  ~  ^^2^^^_^2^-k  ^^ 


X 


2  ^ 

(jO     dT) 


,   X  <  X 


CO 


(19) 


(20) 


_1  X 

v^^)(x)  ~  (u)2_  ^2^^^^  h  ^^^    ^ 


^  /J^A^  d^  + 1 


,  X  >  X^; 


(i) 


v^^\x)   ~  (b^(x)  -  0)2)  ^  exp  A 


lA^ 


,   X  <  X 


To  express  -i^{x,±i^  )    in  terms  of  v^^^  and  v^ '^^  ,  we  note  first  that 


^  Xg  X 


2  ^ 

CO     dT). 


Using  the  relations  between  sinh  and  exponential  functions,  (l5)  may  be 
rewritten  as 


(22)       V  (x,±u)  )~  23  v^3)(^)  _  p-l  ^(M(^ 


where 
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(23) 


P  =  exp 


A  \  /b  (ii)  -  w   dTi  +  7 


For  X  >  X  ,  then,  the  solution  v  can  be  obtained  from  the  connection 
formulas : 


(2k) 


2        2   ~   k 
v^(x,    ±  w)~   (to    -  b    )  2(3   sin 


yoj    -  b    (ti)    dTi  +  ^ 


OJ 


P        cos 


J    yu)   -  b   (t])    dTi  +  I 


To  simplify  the  determination  of  the  Wronskian  we  note  that  from  (ik) 


+  17t 


(25) 


v^(x,±oj)~  e   ^   [±i  v^3)(^)    ^  v^^^(x)] 


,(^), 


Then,    from  (22)       and  (25), 


(26) 


\^(v^,v^)~  e         W 


±iv^3)(^)^^(^)(^)^    2(3v^3)(^)    _    p-l^(^)(^) 


+irt 


W(v^.V2)~   (23±i3"^)e    ^     W(  v^  ^^  3)  ^^^(gpiip"^) 


Tirt 

1+ 


The  final  result  is  obtained  from  (17)  and  (19)  or  from  (18)  and  (20). 

!2_ 
Noting  that  3  is  exponentially  large  in  A,  the  quotient,  y  ,  is  formed 


and  the  exponentially  small  part  neglected.   It  follows  that  for  x  >  x 


0) 


(27) 


W(v^^y^) 


"2  ,  2,  . 
u)  -  b  (x) 


_1   ±in 

k        k 
e 


2iA 


exp 


exp 


iA  \  /w  -  b  (t|)  dT)  +  Y" 


to 

X 


-iaT/o^-  h^i 


Tl)  dT] 


m 


to 


-65- 
Appendix  II 

In  this  section,  we  prove  that  the  integral  given  by  (l),  below, 
is  asymptotically  equal  to  the  integral  given  by  (3).   This  result  allows 
us  to  simplify  equations  (U.U),  (  U.  5)  and  (  U.  11)  . 

We  consider  the  integral, 

(30  z 

(1)       1(A)  =  A  \  f(z)g(Az)exp[iA  \  cp(Ti)dTi]  dz. 

If  g(z)  is  a  continuous  function  with  compact  support  (call  the  support 
£)    ;    I  z|  <  z  for  z  e  i)  )  and  9  and  f  have  continuous  derivatives  in  D  , 
then 


(2)       1(A)  =  K  +  o(|)  , 


where 


00 
(3)       K  =  f(o)  J   g(z)exp[icp(o)z]  dz. 


Proof: 


00       ^  z 

(U)    1(A)  -  K  =  A  \   g(Az)|f(z)exp[iA  \  9(Ti)dTi]  -  f  (o)  exp[  iAcp(  o)  z]Uz. 


We  rewrite  {h)    as 


(5)  1(A)  -  K  =  yA)  +  I^(A), 

where 

00  ^  z  ' 

(6)  12(A)  =  A  \   g(Az)f(z)iexp[iA  I  cp(Ti)dTi]-  exp[  iAcp(o)  z]  Idzj 
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00 

(7)  I   (A)=  aV  g(Az)exp[iAcp(o)z]{f(z)    -  f(o)}  dz. 

Let 

(8)  F,    =  max   |f(z)|  ,    F^  =  max   |f'(z)|  ,   G  =  max   jgCz)]  ,    <!>  =  max   lcp'(z) 

^       z6f)  '^      z€©  z€f)  zeS 

From  (6)    and  (8) 

z 

(9)  llg^^n    <2AF^Gj    I  sin  I  ^    (cp(  t])  -  cp  (  o)  dTi|  dz 

lAz|<z^  ° 

We   introduce  the  variables   ^  and  ^,    defined  by 

(10)  I  =   ATI        ^   =   Az, 
and  rewrite   (9)    as 

c 

(11)  |I2(A)|    <  2  F^G^  |sin  i^(cp(^^)  -  cp(o))d^|dC. 
Since    |  sin  Q\    <   1^|  , 

C 

(12)  1 12(A)  I    IF^G^        I    j  (J^  -  cp(o))d?|dC; 


(13)  l2(^)   <  ^1^^     j  ^    0  -^     d^   dC  =  —^ 

Letting  C  =   Az   in  (7),   we   obtain 

(lU)  13(A)    =y  g(C)exp[i(p(o)z]|f[y   -  f(o)|d^; 
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(15)  |I3(^)|<G^         |ff^j-  f(o)|dC<  G    y  F 


2     A     ^^• 


a<-o  1^1^^ 


Then 

c   "F 

(16)  113(^)1  <    -^-l 

From  (h)  ,    (13)    and  (16), 


(17)  IKA)    -  K|    <     I     , 


so  that 


(18)  1(A)    =  K  +  oJy 
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Appendix  III 

In  this  appendix  we  shall  show  that  the  branch  cut  integral  in 
(9.18)  is  asymptotically  zero.   Also,  the  integral  on  a  semicircle 
in  the  lower  half  of  the  w  plane  tends  to  zero  with  increasing  radius. 

In  order  to  determine  the  branch  points  and  branch  cuts  of  the 

integrand  of  (9- 18)  we  must  examine  the  functions  cp(x,to)  and  \|r(x,cj) 

defined  by  (9.8  )•   For  b(x)  <  co  <  b  ,  ii{x,i^)  >   0  and  arg  i)f(x,u))  =  0. 

When  Lo  is  near  b(x),  by  integrating  by  parts  in  (9.12)  we  find  that 

3/2 

(1)   ^^^'^^^   ^-b(x)V(l)    '      '^-^^-)- 

Hence  \J/  has  a  branch  point  of  order  3/2  at  ib(x).   At  uJ  =  lb  and 

(jj  =  lb  the  functiorP  cp  and  ^   can  be  analytically  continued  by  power 

series  and  therefore  have  no  singularity.   The  functions  ?(x,  1^) , 

cp(x^,w),  q)(x  ,w)  and  their  derivatives  with  respect  to  x  appear  in 
n  O 

(9.16).   We  therefore  have  branch  points  at  ±b(x)  ,  ±b(x  )  and  +b  =  +bfx  ) 

\  /  7    \  q/       o    ^  B 

and  take  the  branch  cuts  to  extend  from  -b(x)  to  b(x),  -b(x  )  to  b(x  ), 

and  -b  to  b  .   The  integral  around  branch  cuts  is  then  an  integral 

enclosing  the  real  oj  axis  from  -b  to  b  . 

00 

From  (1)  we  see  that  \|r(x,uj),  x  =  x,  x  ,  or  x  ,  has  a  branch  point 

O  D 

of  order  3/2  at  x.   Let  us  consider  x  =  x.   For  w  >  b(x ),  \|f(x,Lo)   >  0 
and  :p(x,uj)  >  0;  arg  \|r  =  arg  cp  =  0.   From  (1),  for  b  <  w  <  b(x)  and  above 
the  cut,  arg  ijr  =  ^  and,  from  (9.8  ),  arg  cp  -  -n.      Below  the  cut,  arg  i)a  =  -  ■~- 
and  arg  cp  =  -jr.   Hence  for  -b  <  w<b(x)  arg  '4>   differs  by  2v:   above  and 
below  the  cut.   Since  the  Airy  function  is  entire,  Ai(  X  '  fe   '4)(x,w))   and 
Ai(x  '    e"^   T(x,w))  are  the  same  above  and  below  the  cut.   Since  arg  ? 
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does  not  change  at  b^,  this  must  be  true  also  for  -b  <  co  <  b  and 

for  -b(x)  <  w  <  b(x).   By  a  similar  analysis,  for  -b  <  w  <-b(x),  arg  '^ 

differs  by  kn   above  and  below  the  cut  and  hence  again  the  Aity  functions 

are  the  same.   A  similar  analysis  can  be  carried  out  for  x  =  x  or  x  . 

o     B 

To  analyze  the  derivative  with  respect  to  x,  we  note  first  that  from 
(1)  and  (9.  8) 


(2)     'H^,.):,    (|)2/3  ;  -^  (^)  ,   ;   u.«b(x). 

(-b(x)b'(x)l^/-^ 


By   differentiating  (9. 8  )  with  respect  to  x,  it  follows  that 


(3)  ?'(x,co)[cp(x»]^/2  ^  ^^2_  ^2^^^-jl/2^ 

From  (2)  and  (3)  it  follows  that  arg  ?' (x,uj)  =  0.      This  analysis  also 
holds  for  5c  =  x^,  x^.   Consequently,  since  the  paths  of  integration 
above  and  below  the  cut  are  in  opposite  directions  it  is  clear  that 
the  resulting  integral  is  zero. 

To  analyze  the  integral  on  the  semicircle  in  the  lower  half  plane, 

2   2 
we  note  first  that  for  b  <  w  (real)  there  are  no  turning  points.   On 

the  real  line  then,  the  solution  is  given  by  v(x,  +w)  in  (2.19)  with 

R_i_  replaced  by  the  quotient: 

-   .  ^     2  _ irt        -.  _^ 

(U)  R  =  n   Ai[A"^"?(x^,oj)J  ^  QA±[\^e   ^   9(Xg,co)]   [>   . 
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On  a  semicircle  of  sufficiently  large  radius  in  the  lower  half  w-plane 
it  can  he  shown  that  R   is  bounded.   In  fact,  to  show  that  the 
integral  on  this  semicircle  is  zero  we  need  only  show  that  the  exponents 


X 

X    o 


1/2 


(5)  e|^^  =  -   /+  /  p_b2(Ti)j     dTi  -  ^ 


and 


(6: 


s 


.(2)  _ 


w  -b  (t) 


1/2 


dr)  -  cot 


have  a  positive  imaginary  part  for  Im  w  <  0.   In  this  case,  the  integrand 
of  (9.18)  will  exponentially  decay  in  the  lower  half  plane  as  |w|  ->  ». 
As  a  sample,  let  us  consider 
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X    o 


2  2   V/^ 
oj  -b  (t]  )     dT)-urt, 


(7)  el^'  =  /. 

From  (9.28)  for  the  region  of  interest 


Xg       B 
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We 


note  that  T^-  tg>  0,  so  that   Im[-  w(T   -t    )]  >  0.      The  integral,    I, 
appearing  in   (8)   can  be  rewritten  as 


X  or 


(9)     I  =  oj  /  + 


^B     ^B 


^lll 


1/2 


IV  b^(T 


1/2 


V  ^^(^) 


-1/2 


dJ], 


As      |w|     ^00 


(10) 
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b!(TL 


1/2 


rf   1   - 


2        2 


so  that 


(11)      I 


^B       ^B 


(h2-h2(,)) 


1/2 


IbiLlL 


o     2        2 
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2  ,  2 
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'2 


(b;-b''(Ti))  dTi, 


But 


(12) 
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so  that   for     Im  uj  <  0 
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(13)   Im  I  >  0 

for  |co|  sufficiently  large.   Hence  Im  e_;   >  0  and  Re  ie   <  0  in  the 
lower  tialf  plane.  A  similar  analysis  can  be  carried  out  for  all 
exponents  in  (5)  and  (6).   We  then  conclude  that  asymptotically  the 
integrand  of  (9.  18)  decays  exponentially  as  |w|  ->  oo  in  the  lower  half- 
plane.   On  the  real  axis  the  integrand  also  approaches  zero  as  |lo|  ->  oo. 
With  this  information  it  can  be  shown  that  the  integral  on  a  semi -circle 
in  the  lower  half  plane  approaches  zero  as  the  radius  of  the  circle 
approaches  oo. 
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Fig.l 


The  ray  picture  for  the  case  of  rapidly  varying  initial  data. 

Primary,   tiarned  and  reflected  rays. 

Space-time  diffracted  rays. 
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Fig.    2        b(x) 
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